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Editors’ Foreword 


Experimentation is essential in engineering. But despite the growth of 
elaborate undergraduate laboratories, experimental tuition has rarely kept 
pace with the increasing theoretical content of engineering courses. This 
series of books grew from a belief that the process of learning in engineering 
laboratories would be enhanced by the ready availability of comment and 
advice such as is provided in these texts. Although intended primarily for 
use in universities and colleges of technology, certain sections may also 
serve as more general references to experimental practice. SI units have been 
adopted throughout, supplemented in some instances by the British system of 
units. 

Volume I contains an introduction to the general requirements of experi- 
mental method, and describes the principles and applications of some modern 
instruments. Volume II (Applied Mechanics) and Volume III (Thermo- 
dynamics, Fluid Mechanics, and Heat Transfer) outline sets of experiments 
which range from the simple to the difficult, and from the first to the final 
year of an undergraduate course. The outlines are not, of course, intended to 
be followed slavishly but to offer material on which coherent laboratory 
courses can be based. 

We are indebted to many sources of information, which are recorded in 
the successive volumes. We record also our appreciation of the patience 
shown by the publisher in what soon became a complex task. But most of all 
our thanks are due to the authors for the way in which they have co-operated 
to realize the original idea. It is they who have written the series and any 
virtues it may possess are due to them. 


1971 L.M. 
S.D.P. 
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Preface to Volume [| 


Part I of this volume describes experiments associated with syllabuses in the 
subject normally called ‘Mechanics of Machines’. With reference to a three- 
year full-time honours course, the first eight experiments are intended for the 
first year, the next six for the second year and the last four for the final 
year, though the boundaries and the order of presentation may be changed as 
required. 

It is now quite normal to include an introductory study of vibrations in 
courses on machines and, more recently, automatic control and analogue 
computation have been added. For the first year we have included simple 
experiments to assist in the understanding and engineering applications of 
mechanisms, epicyclic gears, friction, balancing and free vibration, and 
others to illustrate the momentum and energy principles, moments of inertia 
and dynamic loading. For the second year, more ambitious experiments 
take vibrations a stage further, deal with gyroscopic precession and shaft 
whirling and introduce some applications of automatic control. The last four 
experiments assume that the reader has acquired a good understanding 
of principles and a reasonable knowledge of laboratory methods. These 
more-involved experiments on advanced vibrations, control and analogue 
computation are outlined for the reader to develop according to his own 
interests. 

Part II is concerned with stress analysis, though this has been extended 
slightly by including several experiments on structural analysis. The first 
five experiments are intended for the first year, the next four for the second, 
and then three for the final year. The emphasis is on practical techniques, 
leading the student through simple examples of experimental stress analysis. 
The impact of computers has resulted in an increasing concentration upon 
mathematical techniques in undergraduate courses. However much of 
engineering still depends on empirical relationships and experience, and 
this requires a sound appreciation of such techniques of measurement and 
analysis as are described in this section. 

It is hoped that the actual execution of the majority of the experiments 
described will assist the reader in appreciating the subject and will introduce 
him to the wide variety of engineering applications involved. The reader is 
encouraged to bear in mind, at all stages, the measurement and computation 
errors involved and to try to assess the accuracy of the results which he 
obtains. 


1971 C.R.W. 
A.R.L. 
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Simple Mechanisms 





(a) Crank-slider 


Theory 


Kinematic mechanisms with turning and sliding pairs of mechanical elements 
appear commonly in engineering practice. One of particular importance is 
the four-bar linkage with three turning and four sliding pairs in the form of 
the so-called ‘crank-slider mechanism’ illustrated in Figure 1.1. 


Crank A,-~ : Soe , Connecting rod B 








Slider D 


ee Fixed frame C 
7 SS ee 


Ficure 1.1 Crank-slider mechanism. 


The crank A of radius r rotates about a fixed point O in the frame C, O 
being regarded as the origin. Its angular displacement 6 is measured as 
shown in the figure. The connecting rod B of length / is pivoted to the crank 
and to the slider D, the latter being constrained to move along a straight 
line in the frame through O. The motion is two-dimensional in the plane of 
the paper. The displacement of D (relative to O) is given by 


1 1 
= Be ren hee te ith Oa 
x=1tr|cos6 pee 6 893 sin*@—.. | (1.1) 


where gir (1.2) 
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In most practical applications, 6 is assumed to be constant and equal to 
w, where «/2m is the frequency of rotation of the crank. In this case the 
velocity and acceleration of D are given by 


; E tome [ee : 
v= —or| sin 6+ gn 26 + 83 (sin 20 —4sin 46) +. . | (1.3) 


1 
i= —o%r| cos 0+ re 26 + iz (cos 20 — cos 40) +.. | (1.4) 


Unless g is very small, the expansions given in equations (1.1), (1.3) and 
(1.4) converge rapidly and, for g in the common range 3 to 4, it is normal 
practice to re-write (1.3) and (1.4) as 


; : ie 

5 om or sin 6+ 3g sin 26] CS} 
1 

ox —cr4| 00s 0+ F cos 26 (1.6) 


For q very large, the second terms in the above expressions become negligibly 
small, so that the motion of D approximates to simple harmonic. 


Experimental Details 


Given a crank-slider with a typical value of g of about 3, a useful experiment 
would be to investigate to what extent the motion of D differs from simple 
harmonic. Using equation (1.1), the theoretical comparison for displacement 
is shown in Figure 1.2. The greatest difference occurs at 0 = 90° and 270°. 





(@) 30 180 270 360 
@ (degrees) 


FicureE 1,2 Motion of slider compared with simple harmonic motion. 


It would be adequate to measure x for increments in 6 of 30° between 
6 = 0 and 0 =360°. Provided the test mechanism is manufactured reasonably 
accurately, any disagreement with theory should arise solely from inaccurate 
measurement of 6 and x. 
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Discussion 


Experimental verification of equations (1.3) and (1.4) or their approximate 
forms (1.5) and (1.6) would be difficult. Would you, however, accept these 
results having confirmed the validity of equation (1.1)? 


' (b) Hooke’s Joint 
Theory 
Rotary motion is often transmitted between a pair of shafts, the axes of 


which are non-parallel but intersect, by means of a Hooke’s joint, as shown 
diagrammatically in Figure 1.3. 





Ficure 1.3 Hooke’s joint. 


For a pair of shafts out of line by an angle «, the relation between their 
angles of rotation is 


tan¢d = tan 0 cos « (1.7) 


where 6 and ¢ are the angles of rotation of, say, the driver and the driven 
shaft, and are measured from the position where the arm of the ‘cross’ 
associated with the driver lies in the plane of the angle «. The angular speeds 
of the shafts are related by 


; COs a 
oe 1 — sin? 0 sin? « (1.8) 


For constant angular speed 6 of the driver, the acceleration of the driven 
shaft is 


cos « sin? « sin 20 : 
a 1 — sin? 6 sin? « (1.9) 


This is zero only when 0 = 0, 90°, 180°, etc. because when the input speed is 
constant the output speed fluctuates during each revolution. 

For « equal to 20°, equation (1.7) enables a graph of (f — 6) against 
6 to be drawn as in Figure 1.4; the curve is for one complete revolution, that 
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for each half-revolution being identical. The speed ratio $/0 is shown in 
the same figure. Although the maximum difference in position compared 
with a constant-ratio transmission is less than 2°, the speed fluctuates by 





360 


(¢ — 6) degrees 


@ (degrees) 


Ficure 1.4 Displacement difference and velocity ratio for Hooke’s joint. 


nearly +7%. Such a speed fluctuation could result in severe vibration 
and excessive stress variation with resulting fatigue in many engineering 
applications. Can you quote any? 


Experimental Details 


A simple experiment may be arranged to measure (¢ — 6) for various values 
of 6 to compare with the first of the curves of Figure 1.4. Increments in 0 of 
15° between 6 =0 and @ = 180° would give a satisfactory result; it is 
necessary to take a few points in the range 0 = 180° to 360° to check that 
the motion is repeated. It is evident that, since (6 — 9) is likely to remain 
small, both 6 and ¢ must be measured very accurately, i.e. to within about 
1 minute. Any slackness in the bearings of the joint would lead to 
measurement errors comparable with the values of (¢ — 6). 

A drive may be transmitted between two Hooke’s joints in turn, and if 
a is the same for each and the arms of the crosses associated with each end 
of the intermediate shaft are parallel, the velocity ratio will be constant. 
Where this is essential, e.g. for use in a front-wheel-drive vehicle, the two 
joints may be assembled effectively as a single joint. Variations of this are 
the Rzeppa and Birfield joints. To check such a drive, two similar Hooke’s 
joints may be arranged in series and the above experiment repeated. You 
should confirm that a constant-ratio transmission results when the input and 
output shafts of the assembly are either parallel, i.e. the angle « for one joint 
cancels that for the other, or lie at an angle of 2« to one another. 


~ 


SIMPLE MECHANISMS 7 
Discussion 


Do your tests confirm the behaviour of the mechanism as predicted by the 

theory? Do you consider that the errors between the theoretical and experi- 

mental curves are due entirely to lack of experimental precision and not 

because of any shortcomings in the theory? Can you give a theoretical 
; justification for the last result of the final experiment? 


2 ‘ 


Epicyclic Gears 


One important modern application of epicyclic gears is the self-changing 
gearbox fitted to many road and rail vehicles. This utilizes either a train of 
helical spur gears or an arrangement of angled-bevel gears. Whereas ordinary 
gear pairs are easy to understand from the point of view of speed and torque 
ratios, the methods used for calculating these ratios in epicyclic trains are not 
always convincing and experimental study can be valuable. 


Theory of Epicyclic Gear Trains 


Consider the simple gear shown in Figure 2.1. The sun wheel A meshes with 
one or more planet wheels B which in turn mesh with an internally-toothed 


i 


c 


ex 


cap 


Ficure 2.1 Simple epicyclic gear train. 





wheel C. The planets are pivoted on shafts fixed to the carrier or cage D. 
Suppose that V4, Ng, and Nc are the numbers of teeth on wheels 4, B, and C 
respectively; let the speeds likewise be w,4, ws, and wc. Because the pitch 
must be common, it follows that 


No = Na + 2N3 


With the cage fixed, the speed ratio between A and B is inversely propor- 
tional to the ratio of teeth numbers; thus 


w/z = —N,/Na 
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the negative sign being necessary because of the opposite direction of motion 
of the wheels. When the cage moves at speed wp, it follows that 


(W4 — @p)|(@g — Op) = —N,/Na (2.1) 


In using these formulae, a sign convention must be adopted for positive speed 
and used throughout. We may also write 


(@4 — Wp)|(@e — @p) = —No|Na (2.2) 


in which it should be noted that N, does not appear. 

Compounding of two or more epicyclic trains is commonly practised. 
Consider, as a specific example, the gear found in a modern automobile 
transmission which is shown diagrammatically in Figure 2.2. The output is 





FiGurE 2.2 Compound epicyclic gear train. 


always taken from D. The input may be applied at A or at E. The wheels C 
are keyed together and have 30 teeth. Wheel FE has 55 teeth and hence D 
has 115. Wheel A has 48 teeth. Wheels B are carried in the same carrier F 
as are wheels C but are offset as shown so that they can have a practical num- 
ber of teeth. It should be noted that the actual number of teeth of B is 
immaterial. 

Similar formulae to those used earlier may be written 


(wg — )|(@p — Or) = —No/Nz (2.3) 
and (04 — p)\(@p — Op) = +No/Na (2.4) 


In equation (2.4), a positive sign is used because the drive from A to D 
passes through the extra wheels B. The gearbox is used as follows. 


(a) Fis fixed and E is free. The drive is applied at A. Equation (2.4) gives 
w,/@p = 115/48 = 2:39 
which provides first gear. 
(b) Eis fixed and F is free. The drive is again applied at A. Equation (2.3) 
gives 
—wr|(wp — @r) = —115/55 
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Equation (2.4) gives 
(w4 — @;)[(@p — wz) = 115/48 


We have two equations connecting three variables. By eliminating 
@,, we obtain the following ratio Ron 


@,4/@p = (1/48 + 1/55)/(1/115 + 1/55) = 1-45 
which provides second gear. 


(c) A and £ are locked together and both are driven. The whole turns as 
one unit, thus 


o,/o, = 1 
which provides third gear (direct drive). 
(d) Fis fixed and A is free. The drive is applied at E. From (2.3) we obtain 
@,/Wp = —115/55 = —2-09 
which provides reverse gear. 


It may appear that first and reverse ratios are somewhat small for auto- 
mobile use; however, where this gearbox is used, a hydrokinetic torque 
converter provides a further ratio of 2 when required. 

In practice, input torque is applied to A or E by hydraulic clutches and 
E or Fare fixed by means of band brakes. It is important to know how much 
torque must be applied to E or F to prevent rotation when maximum torque 
is being transmitted through the gearbox. By assuming that the torque ratio 
between input and output is the same as the speed ratio, we may find the 
torque required to hold £ or Fin terms of the input torque, both in magnitude 
and direction. 


Experimental Details 


It is assumed that a compound epicyclic gear is available for test. This should 
be equipped, as shown in Figure 2.3, to enable torques to be measured under 
static conditions. Taking the gear described above as an example, condition 
(b) requires a cord to be wrapped round a pulley connected to A and tensioned 
with a body of known mass. Similar cords in connection with E and D are 
restrained by spring balances. Condition (a) is provided for in a similar 
manner. If the gear enables a reversing condition to be provided, i.e. like 
(d) above, the cord associated with D must apply a restraining torque in the 
reverse direction compared with (a) and (b). 

Count the number of teeth on the various wheels and determine the theore- 
tical speed ratio under various conditions in the manner shown in the 
example. To check the ratios experimentally, lock the appropriate parts and 
measure the angle of rotation of the driver which corresponds to one complete 
revolution of the driven shaft. 
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(b) Second gear (d) Reverse gear 


Ficure 2.3 Measurement of torque reaction. 


Using the equipment for static torque determination, find the torque 
ratio between input and output shafts. Note to what extent this differs 
from the speed ratio under each condition. 


Discussion 


Did your tests confirm that the relationships given in equations (2.1) to (2.4) 
are correct? 

When in motion, the transmission efficiency is given by the torque ratio 
divided by the speed ratio. From the results of your tests, are you able to 
judge how efficient the gearbox would be in use? 

Suggest why the transmission of a road vehicle usually consists either 
of a simple gearbox (with alternative ratios) in conjunction with a dry friction 
clutch or of an epicyclic gearbox combined with a hydrokinetic coupling, i.e. 
why you would not expect to find alternative combinations of these units. 


3 ‘ 


Moments of Inertia 


Moments of inertia frequently appear in the equations of angular motion 
of rigid bodies. For a body of regular shape, its moments of inertia can be 
found from measurements of its mass and physical dimensions, but the shape 
of an engineering component is often too complicated for this approach to be 
feasible. The alternative is to measure the values experimentally. Two such 
experiments, both based on simple vibration theory, are outlined below. 
Their application is wide but, to fix ideas, we shall apply them to the particular 
case of a connecting rod of a reciprocating engine. 


Definitions 


The outline in Figure 3.1 represents an irregularly shaped body in relation to 
three mutually perpendicular axes OX, OY, OZ intersecting at some chosen 


% 


TA, 


Ficure 3.1 Mutually perpendicular axes of a body. 


Zi 


origin O. The body is assumed to be made up of particles of which one of mass 
6m at P is typical. The mass moment of inertia of the body about OZ is 
defined as 


I, = Xr26m (3.1) 


where r, is the perpendicular distance of P from OZ and the summation 
includes all the particles in the body. In the limit this becomes 


I, = frz.dm = §(x? + y?) dm (3.2) 
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If we require the moment of inertia about a parallel axis OZ’ (Figure 3.2), it 
can be shown that J,, and J, are related simply by 


I, = I, + ma (3.3) 


where m is the total mass of the body and a is the perpendicular distance 
_of O'Z' from OZ. Equation (3.3) must be used with care; it applies to all 

shifted but parallel axes O'Z’ only if the centre of mass G of the body lies in 
OZ. 





FiGureE 3,2 Parallel axes. 


Frequently we are concerned only with axes perpendicular to one given 
plane (OXY in our case), and it is convenient to re-write (3.3) as 


I = Ig + ma? (3.4) 


where J is the moment about some axis distant a from the centre of mass and 
Ig is the moment about G itself. 


Theory of the Compound Pendulum 


Consider a rigid body, e.g. the connecting rod shown in Figure 3.3, free to 
swing in a vertical plane about a horizontal pivot at A. For a small angle 
of swing, the free oscillation is governed by the local gravitational field 
according to the equation of motion 


I§ + mga = 0 (3.5) 


where J is the moment of inertia about A, 8 is the angle of swing, m is the 
body mass and a is the distance AG. The periodic time T, of the motion is 
given by 


T, = 2n(I/mga)* (3.6) 
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Note that, since J, m and a are invariable properties of the body, the periodic 
time depends on the value of g at the place where the experiment is performed. 
By equation (3.4), we obtain 


Tq = (mgaT2/4n?) — ma? (3.7) 





FicurE 3.3 Free vibration of a connecting rod. 


Similarly, for small free oscillations about B, we obtain 
Ig = (mgbT2/4x7) — mb? (3.8) 


Since /=a+b, eliminating I, between equations (3.7) and (3.8) we can 
obtain a (or db) in terms of J, i.e. 

_ _— alts — 47°F 

~ gT2 + gT?2 — 87°] Co 


which enables the position of G to be found. (It is assumed, of course, that 
G lies in AB.) 


Theory of the Trifilar Suspension 


Figure 3.4 shows a uniform and circular table suspended in a horizontal 

plane by three equidistant flexible vertical wires of equal length /, each 

being attached to the table at radius r from its centre O. We are interested 

in the free angular vibration of the table about a vertical axis through O. 
For small amplitude oscillation, the equation of motion is 


Ipp + (mgr = 0 (3.10) 


where ¢ is the angular displacement about O and J, and Mr refer to the table. 
The period of the motion is given by 


Tr = 2n(Izl/mygr?)* (3.11) 
and hence Ip = T7mzgr?/477] (3.12) 
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_ Suppose that we now place the connecting rod on the table, with G coin- 
ciding with O, so that it moves with the table. The new period T will be 
related to the resultant moment of inertia J by 


I = T?mgr?/4n?1 (3.13) 
_ Subtracting (3.12) from (3.13), we obtain J, for the connecting rod: 
Iq = (T?2m — T2mz)gr?/4r71 (3.14) 





Ficure 3.4 Trifilar suspension. 


Experimental Details 


For the compound pendulum experiment, equation (3.7) shows that we 
require. values for m, g, a and T,. We can determine m by using a weighing 
apparatus and assume the proper value for g. The reader should devise at 
least two methods of locating G directly in order to determine a. Find T, 
by measuring the time taken for a number of cycles after having set the 
connecting rod swinging on a knife-edge at A with small amplitude. This 
should be done several times, for different numbers of cycles, to determine the 
range of experimental variation in T,. After selecting the best value for Ty, 
I, can be calculated from (3.7). 

Alternatively T, can be found after setting the connecting rod swinging 
about B. Then, after measuring /, a can be found from (3.9) and inserted in 
(3.7). This determination of a should confirm the earlier direct measurement 
of a. 

The trifilar suspension test consists of determining the periodic times, in a 
manner similar to the above, with and without the connecting rod in position 
on the table. Measure the mean length / of the three wires and check that 
this does not alter significantly with the added weight of the connecting rod. 
Also measure the mean value of r. Remove the wires from the table and mea- 
sure my with a weighing apparatus. In determining J, from equation (3.14), 
remember that m is the combined mass of table and connecting rod. 
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Discussion 


State the major assumptions in the theory leading to equations (3.6) and 
(3.11) and discuss how well these were satisfied in your experiments. 

Comment on the main sources of error in each experiment and estimate 
their effect on the J, values found, i.e. estimate the percentage accuracy in 
each value. 

What are the important factors in the design of the experimental apparatus ? 
For example, must (a) the knife edge for the pendulum be sharp and truly 
horizontal, (b) the trifilar wires be truly flexible and of negligible mass, (c) 
1> r in the trifilar suspension? 
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Friction in Machines 





Friction is both the enemy and the friend of the engineer. In this composite 
experiment we shall investigate the frictional resistance caused by conven- 
tional bearings, a source of power loss in rotating machinery. We shall 
also study the deliberate use of friction in clutches and belt drives. 


Theoretical Assumptions 


Consider a body of moment of inertia J rotating freely at angular speed w 
and apply a resisting torque T. The equation for the subsequent motion is 


Ido|dt = —Tp (4.1) 


It is normal to assume that 7, takes one of two forms. ‘Dry’ or Coulomb 
friction assumes that 7’, is constant, and it follows that the speed of the 
body decreases at a constant rate. Viscous friction occurs when T, varies 
directly with w. Let Tz = aw where a is a constant. Substitution in (4.1) 
and integrating gives 


w = w, exp {—(a/D)t} (4.2) 


which means that w decreases exponentially with time from the initial 
speed @. 


Belt Theory 


When power is transmitted from one shaft to another by means of a dry 
belt in tension, it is assumed that Coulomb friction acts between the contact 
surfaces of the belt and pulleys. This gives rise to the idea of a coefficient of 
friction relating the maximum possible friction force to the normal force. 
When motion is transmitted in the normal way, the ratio of these forces is 
less than and we assume that no slipping occurs between the surfaces. 
A knowledge of y enables us to estimate the maximum power which can be 
transmitted. If we prevent the driven shaft from rotating, the belt slips on 


one or both of the pulleys, and under these circumstances we assume that the 
ratio of the forces is equal to yu. 
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In Figure 4.1, the tensions at each end of the stationary flat belt are F, 
and F», and the belt is in contact with the pulley over an angle 0, the angle of 
lap. The pulley, of radius r, is rotating at constant speed w. Simple flat-belt 
theory relates the tensions by the following 


F,/F. = exp “0 F (4.3) 


If the inertia associated with the pulley shaft is 7, equation (4.1) applies. 
The resisting torque afforded by the belt is given by 


Tr = (Fh, — F)r (4.4) 
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Ficure 4.1 Friction between a belt and a pulley. 


We assume that F, and F, (and therefore Tz) remain constant as w decreases. 
Eliminating F. from (4.3) and (4.4) and substituting in (4.1) gives 


Ido/dt = —F,[l — exp (— p6)]r 
Integrating, I(@, — @2) = Fyr[l — exp ( — 6)](t2 — 41) (4.5) 


where , and wz, are the speeds at times ¢, and f,. Suppose the initial speed is 
; hence t; = 0. Let the speed fall to zero at t,; hence w, = 0. Rearrange- 
ment of (4.5) enables us to write 


# = [—In(1 — Jo,/rF,t2)}/6 (4.6) 


This shows that for constant 4 the expression in the square brackets varies 
directly with 6. Note that F, does not appear explicitly. 
For a V-shaped belt, of included angle 2«, (4.3) is modified to 


F,/F. = exp 6 cosec « (4.7) 
and hence (4.6) becomes 


# = [{-ln (1 — Ja,/rF,t.)} sin «)/6 (4.8) 


Clutch Theory 


In Figure 4.2, torque is transmitted between two shafts in line by using a 
plate clutch in which contact occurs between radii r, and rz. Two distinct 
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theories exist. The first assumes that the pressure is uniform over the area of 
contact; this results in the relation 


Ty = 3uP(ri — Qi — 12) (4.9) 


where P is the axial force between the plates and yu is the friction coefficient 
_ acting when slip occurs under torque Tp. 


& 





Ficure 4.2 Friction in a plate clutch. 


The second theory assumes ‘uniform wear’ which is taken to mean that 
the pressure varies inversely with the radius. This results in the alternative 
relation 


Tp = 3¢P(y + 12) (4.10) 


Equation (4.9) gives the larger value for Tz but the difference is small except 
when r, <p. 


Experimental Details 


We shall envisage a composite test rig on the lines of Figure 4.3. A d.c. 
electric motor A, having a maximum power output of about 1 kW at 2000 





Ficure 4.3 Apparatus for studying friction in machines. 


rev/min, drives a gearbox B of speed ratio 0-25. The output shaft drives one 
side of a single-plate clutch C. Axial force is applied to the clutch by a spring 
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balance operating through an equal-arm bell-crank lever D. The clutch 
may be engaged by hand by turning a second bell-crank lever E. 

The clutch output shaft carries a massive body F and two pulleys, one 
flat and one V-shaped. The appropriate belt may be wrapped round either 
pulley by apparatus at G, a spring balance being provided to measure the 
larger tension F,; the angle of lap 6 may be varied in discrete steps of 30° 
from zero to 180°. The speed of the output shaft may be measured by a hand 
tachometer applied to the end of the shaft. 

(a) Bearing Friction and Windage. After starting the motor, the output 
shaft is run up to full speed by hand operation of the clutch with both spring 
balances slackened off. The clutch is then disengaged and the times for the 
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Ficure 4.4 Bearing friction and windage. 


output shaft speed to fall successively to 400, 300, 200, 100 and zero rev/min 
are found with the aid of the tachometer and a stop watch. This should be 
repeated to check consistency. A curve is drawn through the readings plotted, 
as in Figure 4.4. 

To compare with equation (4.1) (for T, constant) and with (4.2), a curve 
for each law is shown plotted in Figure 4.4. In drawing these curves the 
time taken for the speed to fall to half of the initial speed is assumed to be 
common. It is evident that T, in the test depends on both Coulomb and 
viscous friction. 

(b) Belt Investigation. Starting with the flat belt, for various combinations 
of 6 and F, (as set on the spring balance) the output shaft is run up to full 
speed by using the clutch. The speed w, is read and the clutch disengaged: 


«~ 
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the time ft, taken by the shaft to stop rotating is measured. It should be ob- 
served that F, remains sensibly constant during the deceleration and that 
the deceleration is constant. 

Given that J for the output shaft is 0-504 kg m? and that the mean diameter 
of the flat belt is 18-1 cm, the numerator of the right-hand side of (4.6) can 
be calculated for each run. Some typical results are tabulated below. 





6 (rad) 7/6 = af8 [8 p22 ae/3 5S x/6 Sx/6 0 ot 
F, (N) 100 100 40 60 80 40 80 40 
@, (rev/min) 498 486 495 496 492 493 483 492 
ta (S) ht AON TKS 60 34> 601 30.253 
—In (i — I@,/rF,t.) 0:18 0:34 0:37 O51 0-75 0-92 0-89 1-13 


These results are plotted in Figure 4.5. It is reasonable to draw a straight 
line through the test points (to intersect the origin), the slope of which is the 
coefficient of w. 

The test is repeated using the V-belt and pulley. It is given that the belt 
mean diameter is 16-8cm and « = 19°. A similar set of results gives the 
second curve of Figure 4.5, from which » for the V-belt can be measured. 

Note that it has been assumed that bearing friction plays a negligible 
part in decelerating the shaft during the belt tests. Comparison of the times 
taken with those in test (a) shows this assumption to be reasonable. 

(c) Clutch Investigation. In this test, a stroboscope is used to illuminate the 
edges of the clutch plates which are marked with a pattern in contrasting 
colours. When the flash frequency is synchronized with the shaft speed, it is 
easy to tell whether or not the clutch is slipping. In this test the flat belt is used 
as a brake. Choosing 6 = 47 and using the u-value found in the previous 
test, the proportionality between braking torque T, and F, on the spring 
balance can be established as 


Tr = kF, (4.11) 


The motor is started and the axial load P on the clutch is set between 
80 N and 280 N by means of the spring balance at the clutch. F, is increased 
until the clutch is observed to slip. A typical curve of F, against P is shown in 
Figure 4.6. This is evidently a straight line, the slope of which relates F, to P. 
Using (4.11) and (4.9) and given that the external and internal diameters of the 
clutch face are 15-9 cm and 10-8 cm respectively, « for the uniform pressure 
theory can be found. Using (4-10) instead of (4-9), 4 can be again found, this 
time for the ‘uniform wear’ theory. 
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Ficure 4.5 Belt friction. 
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FiGurE 4.6 Clutch friction. 
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Discussion 


Check that you can obtain equations (4.2), (4.3), (4.7), (4.9) and (4.10) 
and note the assumptions made. 

Does your bearing friction test confirm Figure 4.4? Comment on the 
accuracy of calculating In (1 — Jw,/rF,t,) for low and high values of 6. Were 
jyou justified in ignoring the bearing friction during the belt tests? 

Confirm that you arrive at the same y-value (to two significant figures) 
for the clutch using both theories. Why are they so similar in this case? Can 
you conclude that either of the theories is the more correct ? Can you explain 
the reason for the intercept on the vertical axis of the curve in Figure 4.6? 


2—II 


5 ‘ 
Balancing of Rotating Bodies 


Lack of balance in machines causes vibration. For various reasons such 
vibration should be avoided or, at least, minimized. This experiment is 
concerned with the balancing of a set of bodies rotating about a common axis. 


Theory 


Consider a small body A of mass m, attached to a rigid shaft by means of a 
light rigid rod so that, when the shaft is constrained by frictionless bearings, 
A rotates at radius r4 as in Figure 5.1. When the shaft is other than vertical, 





Al Le 
mOe 


Ficure 5.1 Static balance. 


because of gravity it will turn until the body A is at its lowest possible posi- 
tion. Clearly it would not turn if another body B was placed opposite A, 
provided that 


Maly = Mr, (5.1) 


It is not necessary for A and B to be in the same plane of rotation. Under 
these circumstances, we say that A and B are in static balance. 

When we rotate the shaft at angular speed w, A exerts a force myr,w? 
radially outwards on the shaft. We can again balance A by means of a body B 
according to equation (5.1), since w is common. However, it is now essential 
for A and B to be in the same plane. This is elementary dynamic balance. Note 
that static and dynamic balance depend on gravitational and centrifugal effects, 
respectively. 


For A and B situated in planes / apart, the equal dynamic forces give rise 
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to a couple m,r,/w? in a plane containing the shaft and rotating therewith. 
Complete dynamic balance is not obtained unless an equal and opposite 
couple is provided. 

We must now replace B by two bodies C and D opposite to A, as shown in 
Figure 5.2. For the centrifugal forces to balance, since w? is common, 


Mar, = Melo + Mylp (5.2) 
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Figure 5.2 Dynamic balance. 


za 


Now take moments in the plane of the bodies about any point O in the shaft; 
for dynamic balance the resultant moment must vanish. Since w* is common, 
Mal ala — Mcrcle — Myf plp =z.0 (63) 


Obviously we can place O where we like in the shaft but, should it then be 
necessary to measure one of the /’s in a direction opposite to the others, its 
sign must be changed in equation (5.3). It is often convenient to eliminate one 





Ficure 5.3 General situation for dynamic balance. 


of the terms in equations like (5.3) by placing O in the plane of the appropriate 
body. Both equations (5.2) and (5.3) must be satisfied for dynamic balance; it 
follows that static balance is achieved at the same time. 

To study the dynamic balance of a large number of rotating bodies, 
we require a more general theory leading to the construction of force and 
moment polygons. The bodies need not be set opposite one another, as in 
Figures 5.1 and 5.2, but may be arranged as shown in Figure 5.3. Consider 
bodies A, B, C,. . ., N, set at angles 6,, 92, etc. from some datum direction, 
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rotating in planes distant /4, /,, etc. along the shaft from some reference plane 
at O. 

It is now convenient to represent the centrifugal forces by vectors and to 
draw a force polygon with component vectors proportional to mr in magnitude. 
Instead of equation (5.2), dynamic balance requires that we obtain a closed 
polygon, as in Figure 5.4. It does not matter whether we draw each vector in 
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FiGuRE 5.4 Force and moment polygons. 


direction radially outwards or radially inwards so long as we follow the 
same rule throughout. 

The moments can likewise be represented by vectors. Now the vector repre- 
sentation of a moment is a line drawn normally to the plane of the moment. 
Such vectors give rise to a moment polygon, as in Figure 5.4. Strictly, each 
moment vector should be perpendicular to the corresponding force vector, 
but there is no need to make this distinction. To correspond to equation (5.3) 
it is only necessary for the moment polygon to close. Where the /’s are mea- 
sured on both sides of the reference plane, the direction of the moment vec- 
tors must be reversed accordingly. It is usually convenient to take the reference 
plane at one of the bodies for which m, r and/or 6 is unknown; this eliminates 
the corresponding vector from the moment polygon. 


Experimental Details 


A suitable apparatus consists of a horizontal shaft carried in bearings attached 
to a frame which is suspended by light springs, as shown in Figure 5.5. 
Keyed to the shaft are four (or more) discs, the whole being in dynamic 
balance. The shaft is driven by an electric motor via a flexible drive at an 
angular speed @ below its whirling speed but such that w/27 is above the 
resonant frequency of the frame on its suspension (see later experiments). 

The discs are slotted to permit the fixing of bodies of various mass at 
various r and 6 by means of bolts and knurled nuts. Any lack of dynamic 
balance introduced is shown by movement of the frame on its suspension. 

First, to investigate the static balance, remove the flexible drive and fix a 
body A to one of the discs. Note that the shaft rotates under gravity until 
A is in its lowest position. Now fix a body B opposite to A, both to the same 
and to another disc, such that equation (5.1) is satisfied. Note that the shaft 
will remain in equilibrium in any position. 
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With B on the same disc as A, re-connect the flexible drive and switch 
on the motor. Note dynamic balance which is shown by the lack of movement 
of the frame on its suspension. 

Using three different discs, fix two bodies C and D both opposite to A, 
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FiGure 5.5 Dynamic balance machine. 


at radii which obey both equations (5.2) and (5.3). Choose values for m andr 
which permit C and D to be (a) one either side of A, (b) both on the same 
side of A. In both cases switch on the motor and check the dynamic balance. 

Now suppose that you are given the information in the following table. 








Plane m (kg) r (cm) 60 1(cm) mr mrl 
A 0 
B 1:0 4:0 180° 7:4 
€ 265° 13-4 
D 40° 20-8 





Note that it is assumed that A, B, C, D are placed in order along the shaft and 
that the plane of A is regarded as the reference plane. You are required to find 
values for m4, Mc, Mp, a, Tc, Vp and 0,4. Each vector polygon will reveal 
two unknowns, leaving some arbitrary choice for the remaining quantities. 

Since mrl for B is known, first draw the closed moment polygon, using 
the given angles; because mr/ for A is zero, you do not need 6). This will 
reveal mri for C and D. Now calculate mr for B, C and D and draw the closed 
force polygon. This will reveal mr and 6 for 4. You now have some choice in 


the selection of the unknown m and r values. 
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Fix bodies to the discs according to your completed table and check 
the dynamic balance. 

As a further example, complete the following table with the aid of vector 
diagrams and confirm your results by a test. 


S 
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Plane m (kg) r (cm) 0 1 (cm) mr mel 
A 10° 0 
B 1:5 3-0 TA 
C 1-0 4-0 13-4 
D 1:0 3-0 20:8 
Discussion 


What principles of mechanics are invoked in the theory of dynamic balanc- 
ing? Why is it possible to draw the vectors of the moment polygon in the 
same directions as those of the force polygon? 

Note any shortcomings in the experimental equipment. Did you take 
into account the different lateral displacements of the centres of mass of 
different sized bodies relative to their respective discs? Does it matter that 
the bodies used in the tests are of finite size? 


Note on the Balancing of Reciprocating Bodies 


In an earlier experiment we considered the crank-slider mechanism. From 
equation (1.6), for a reciprocating body of mass mp, the accelerating force 
at the body is very nearly —m,w?r{cos 6 + (cos 20)/q}. This consists of 
two forces, the first, —m,w?r cos 6 being called the primary force, and the 





Ficure 5.6 Primary and secondary cranks. 


second, (—m,w?r cos 26)/q being called the secondary force. The reactions of 
these forces on the crankshaft are equal and opposite. 

We may now visualize a body P rotating with the crank and another Q 
also rotating about O but at twice the speed of the crank, as depicted in 
Figure 5.6. These could be arranged to provide exactly the same variation of 
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reaction on the crankshaft in the direction OD, as is caused by the motion of 
the reciprocating body considered above. This is why problems on reciprocat- 
ing balance are similar to those on rotating balance except that two sets of 
polygons must be drawn. By equation (1.4), forces exist in addition to the 
primary and secondary forces, but these are almost always too small to matter. 

In a reciprocating engine, mz is taken as the mass of the piston plus some 
portion of the mass of the connecting rod. The remaining connecting rod 
mass is taken as acting at the crank pin, where it can be regarded as a simple 
rotating body. 


6 ‘ 


Free Vibration of a 
Simple System 


The natural frequency or frequencies of vibrating systems are of great impor- 
tance and must often be determined experimentally. A system vibrates in 
its natural or free mode(s) when any disturbing agents from without the 
system are zero or constant in magnitude. A familiar example is the vertical 
motion of a body suspended by a tension spring from a fixed support; it is 
necessary to disturb the system in order to set the body in motion, but we are 
interested in the free motion after the disturbance has been removed. 

The example chosen for the present experiment is a one-degree-of-freedom 
linear torsional system. 


Theory 


Consider the system shown in Figure 6.1, which comprises a body suspended 
at its centre of mass by a light rod or wire from a fixed support. This system 
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Ficure 6.1 Simple torsional system. 
can vibrate freely in a number of ways. We could cause the body (a) to 
swing from side to side like a pendulum, or (b) to move bodily in a circular 
horizontal path without rotating, or (c) to rotate about its centre of mass 
causing the wire totwist.We say that the system has several degrees of freedom. 


We are interested only in (c) above and we must restrain the system from 
vibrating in any other manner, 
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Suppose that the wire is circular and of uniform diameter d and length /. 
From the simple theory of torsion (see Experiment 20), 


T|J = Goll (6.1) 


where @ is the angle of rotation of the body and is therefore the angle of 

* twist over the whole length of the wire, T is the torque imposed on the 
wire by the body (and by the support), G is the modulus of rigidity of the 
material of the wire, and J is the second moment of area of the wire about 
the axis of twist and is equal to 7d*/32. We are interested in the stiffness k of 
the wire, i.e. the torque per unit angle of twist; hence 


k = 7/0 = GJ|l = Grd*/321 (6.2) 


For convenience, take the body as being a circular disc of mass moment 
of inertia J about its point of support (centre of mass). The reader should note 
the essential difference between J and J. Suppose that a cylinder fitted under 
the body rotates within a fixed cylinder, the space between the cylinders being 
filled with a viscous medium, Figure 6.2. This constitutes a dashpot, the 





FicureE 6.2 Damped torsional system. 


viscous friction coefficient of which we suppose is ¢, i.e. c is the resisting torque 
on the body per unit angular speed. Then, equating torques for some arbitrary 
6, we obtain the system equation of motion 


16+ch6+k6=0 (6.3) 


This is a second-order homogeneous linear differential equation in 0. 
It has several solutions depending on the relative magnitudes of J, c and k. 
For c2 < 4/k, when the body is released from an angle 6 measured from the 
equilibrium position, we obtain the damped oscillation 


6 = O(1 + a/b?) e* cos (bt + B) (6.4) 
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where a is the damping factor, b/27 the natural damped frequency and 
the damping angle, given by 


a= —¢/2I 
b = (k/I — c2/412)t —. (6.5) 
B = arctan (a/b) 


Two other cases are of interest. For no damping, c = 0 and instead of 
(6.4), we obtain the undamped oscillation 


6 = 6 cos b,t (6.6) 
where the natural undamped frequency, b,/27 is given by, 
b, = (k/D (6.7) 
For c? = 4Jk, we obtain the aperiodic solution 
6 = Ge (1 — at) (6.8) 


The three results are compared in Figure 6.3. For c = 0, the motion is 
said to be undamped; for c? < 4Jk, it is underdamped; for c? = 4Ik, it is 





FicureE 6.3 Free motion of the system. 


critically-damped. There is a fourth case known as overdamped, when 
c? > 4/k, but this is of no interest in the present experiment. 


Experimental Details 


This follows the lines of Figure 6.2. In a typical case, the suspension wire is of 
steel, 0-90 mm in diameter and 0-91 m long. The modulus of rigidity is 
81:5 x 10° N/m?. The disc, with cylinder attached, is of brass and is well- 
balanced about its point of attachment to the wire. Its moment of inertia is 
0-002 64 kg m? and the cylinder diameter is about 10 cm. The upper surface 
of the disc is graduated in degrees around the periphery. A means is provided 
for releasing the disc from an angular position displaced from the equilibrium 
position; for convenience, # is made 180°, and it is essential that this does not 
twist the wire beyond its elastic limit. 
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The fixed cylinder is arranged so that it can be moved vertically relative to 
the moving cylinder. The radial clearance must remain uniform and equal to 
about 0-5 cm. This cylinder is filled with a mixture of fuel oil and lubricating 
oil gauged so that, over the range of vertical movement of the cylinder, the 
effective viscous damping coefficient covers most of the range 


0< c? < 4lk 


i.e. so that the free motion of the system can be investigated over most of the 
underdamped region. 

Practise releasing the disc from its displaced motion until you are able to 
set it in smooth angular motion. Check that the equilibrium position of the 
disc corresponds to 0° on the scale. 
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FiGureE 6.4 Free underdamped response. 


Lower the fixed cylinder until the moving cylinder is clear of the oil. Start 
the vibration and record the amplitude each time the disc comes to rest, 
until about ten complete cycles have been completed. This can be done easily 
and accurately after a little practice, especially if one person reads whilst 
another records. Repeat the experiment to check consistency and measure, 
with a stop watch, the total time for the ten cycles. Find 5, and compare with 
the value calculated from equations (6.7) and (6.2); for the latter, you will 
need to take given values for J and G and to ensure that you have a very 
accurate value of d. It is to be expected that the free motion will be slightly 
damped even though no oil damping is used. 

Next adjust the lower cylinder until ¢ is large enough to ensure that the 
amplitude decreases from 180° to about 20° after ten complete cycles. Repeat 
the test as above and estimate b by time-averaging over the whole ten cycles. 
Do not time from the instant of release but from the first occasion when 
the disc stops. Now try to measure b separately during the early and late 
part of the transient; there should be no difference if the system is truly 
linear. 
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Plot the amplitude against peak number (which is proportional to elapsed 
time) following the style of Figure 6.4. Draw smooth curves to represent the 
envelope of the vibration, and sketch in the disc displacement between peaks; 
from equation (6.4), the envelope should be the curves + e*. Also plot the 
logarithm (to base e) of the amplitude against peak number. The envelope 
should now consist of straight lines of slope +a. 

Using equation (6.5), find ¢ from the slope of the logarithmic envelope. 
Also check that the following relation is obeyed: 


b2 = b? + a? (6.9) 


Repeat the above but with about twice the amount of effective damping. 


Discussion 


Did your calculated and test values of b, agree closely? If not, what do you 
consider to be the primary source of error? What caused the slight damping 
which occurred when the dashpot was rendered inoperative? Did this cause 
any significant error in the determination of b,? 

What part does gravity play in this experiment? Does it affect any of the 
results? Why is it necessary to avoid timing between the moment the disc 
is released and when it first comes to rest? : 

Can you assess whether the oil damping was truly viscous or otherwise? 
Can you obtain equation (6.9) from the theory given earlier? Did 6 differ 
greatly from b,? 
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Dynamic Loading of a System 
with One Degree of Freedom 


When a system is in equilibrium with a load imposed upon it, we speak of 
static loading. When potential and/or kinetic energy is transferred to the 
system during the imposition of the load, we have dynamic loading. It is 
important to appreciate that the stresses arising in the system during dynamic 
loading may be many times greater than those occurring under static loading 
for the same load. 


Theory 


Consider the system with one degree of freedom shown in Figure 7.1, which 
consists of a body of mass m and a massless spring of linear stiffness k. Motion ° 


Body of 
mass m 
mg 


h 







Spring of 
stiffness k 


Ficure 7.1 Simple spring-body system. 


is restrained to the vertical direction along the axis of the spring. The lower 
end of the spring is attached to a fixed support and the upper end is free to 
move. Suppose that the deflection of the spring (displacement of the upper 
end) under the action of the load is x, measured as positive in the downward 
direction. 

First suppose that the body is gently lowered until its weight mg is supported 
by the spring in equilibrium. The static deflection is clearly given by 


x, = mglk (7.1) 
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The stress in the spring is proportional to x, say equal to px; hence the static 
spring stress is of magnitude pmg/k. 

Now suppose the body is released to fall freely under gravity through a 
vertical distance h until it reaches the upper end of the spring, and suppose 
that it then attaches itself to the spring. The system will first come to rest 
when the spring has moved through its dynamic deflection xz. Assuming no 
loss of energy, we may equate the strain energy of the spring to the reduction 
of potential energy of the body. 


$kx2 = mg(h + Xa) (7.2) 
Solving for xq gives 
=e ae) | 
x= El 14 (1455 (7.3) 
Consider the special case when h = 0. Equation (7.3) then shows that 
Xq=2mg/k or O (7.4) 


The first value gives the dynamic deflection when the system first comes to 
rest; we note that this is twice the value of x, given by equation (7.1). The 
spring stress is likewise twice the static stress, i.e. 26mg/k. In this position the 
spring force is 2mg, sufficient to return the body towards its initial position. 
In fact, the second value of xq given by equation (7.4) shows that the system 
next comes to rest when the body has returned precisely to its initial position. 

It is now evident that an undamped free oscillation of the system has been 
initiated, the (downward) displacement of the body being 


x = 3x,(1 — cos D,ft) (7.5) 
where, according to equation (6.7), b, is given by 
b, = (k/m)! (7.6) 


The amplitude is x,/2, in this case x,. 

It is instructive to apply the impulse-momentum equation (see Experi- 
ment 8) to this phenomenon. The net downward force on the body F may 
be written 


F= mg — kx 
2. 
= mg —5 — cos b,ft) 
= mg cos D,t (7.2) 


The time integral of F equals the change of momentum of the body, i.e. 
A(mx) = Jmg cos b,t dt (7.8) 


The momentum of the body is zero initially. The body reaches its static 
equilibrium position after one quarter-cycle, ie. when t = 7/2b,. At this 
instant, the momentum is given by equation (7.8) as mg(m/k)*. This shows 
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that the (downward) velocity is g(m/k)*, which may be confirmed by differen- 
tiating equation (7.5). Similarly, after one half-cycle the momentum of the body 
is found to be again zero, which accords with the facts. 

In the general case, when h> 0, the body first comes to rest when 
Xq > 2mg/k, see equation (7.3), the spring stress then being greater than 
when h = 0. It is evident that for large h the stress can be many times greater 
than that occurring under static loading. The subsequent motion of the 
body is again undamped oscillatory motion about the static equilibrium 
position. Its amplitude is now greater than x,, so that during part of each cycle 
the spring is in tension, reversal of stress taking place when the spring 
deflection is negative (i.e. upwards). 

An alternative analysis can be applied when h > 0. When the body has 
fallen through A, its reduction in potential energy given by mgh equals the 
kinetic energy acquired given by 4mx2, assuming that no energy has been 
lost. Instead of equation (7.2), we may write 


$kx2 = 4mx2 + mgxa (7.9) 
Solving for xq gives 
_ me #24)" 
w= [ra (14 24 a. 


Since, in this case, %) = (2gh)?, equation (7.10) is identical to equation (7.3), 

In practice, due to inherent damping effects which have been neglected. 
the amplitude of the motion will gradually decay to zero as energy is trans- 
ferred as heat, leaving the system in its static equilibrium configuration. We 
have, of course, neglected any loss of energy due to impact where this applies. 


Experimental Details 


It is convenient to use a light elastic cantilever rather than a coil spring; 
this is supported in a rigid frame, as shown in Figure 7.2, and the load is 







Cantilever 


Strain gauges 


Ficure 7.2 System for dynamic loading. 
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applied at the free end. In conjunction with a body of mass 50 g, a suitable 
cantilever is one made of high-tensile steel having the following dimensions: 
length 30 cm, width 2 cm, thickness 1 mm. For these parameters, k is about 
38 N/m and x, is about 1-3 cm; the natural period of the system is about 
0:22 s. To avoid overstressing the material of the cantilever, x, should be 
restricted to about 8 cm, which requires h to be limited to about 15.cm. 
Strain gauges (see Experiment 25) are fixed to the upper and lower faces 
of the cantilever near the fixed end; these feed recording apparatus which is 
calibrated to indicate the stress in the cantilever and also the deflection under 
the load (which are proportional to each other). A signal displayed on a 
cathode-ray oscilloscope (C.R.O.) provides a deflection-time record. 


Body 









Cantilever 
Trigger lever 
Split lever 
277 Body 
a Li 3 Cantilever 


FiGure 7.3 Detail of release mechanism. 


The mechanism for releasing the body is shown in detail in Figure 7.3; 
in the plan view, the cantilever is shown displaced from the body for clarity. 
The body is shaped with two ears which are supported by the ends of a 
split lever. The latter is retained against the tension of a spring by the trigger 
lever. When the split lever is released, the ends supporting the body move 
with an acceleration well in excess of g, thus permitting the body to begin 
its fall freely under gravity. 

The underside of the body is fitted with a cork insert which, on contact 
with the cantilever, is penetrated by a pair of sharp spikes attached to the 
cantilever. The friction between the spikes and the cork assists the body and 
cantilever to move together during the subsequent motion. The whole 
release assembly may be moved vertically relative to the frame and fixed 
according to a scale which indicates the value of h. 
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Move the cantilever by hand and check the calibration of the deflection 
indicator and the C.R.O. Fix the body to the cantilever and measure the 
static deflection x,. Cock the release mechanism and offer this up to the body 
until the ends of the split lever just touch the ears; check that the scale reads 
x, (negatively). The scale may then be relied upon to indicate h. Check 
that the spikes on the cantilever penetrate the cork insert in the body 
sufficiently to ensure that the body rests firmly on the cantilever following 
release of the load from any height between zero and 10 cm. 

Start the tests by releasing the load with h set at zero. Note x4 and the stress 
when the system first comes to rest, and observe the subsequent motion. 
Repeat several times to check consistency. Carry out similar tests at about 
eight different values of h between zero and 10 cm. 


Discussion 


Compare the behaviour of the system with that predicted by the theory. 
To what extent do the measured and calculated values of x, differ as h is 
altered? Give reasons for the differences, taking into account 


(a) the energy lost in impact, 

(b) friction and damping effects, 

(c) the mass of the cantilever, and 

(d) the shortcomings of the test equipment, including dynamic attenuation 
of the stress signal. 


8 . 


A Momentum and Energy 
Experiment 





As a consequence of Newton’s Second Law we have the Principle of the 
Conservation of Momentum, which states that the momentum of a system 
of bodies is constant provided that no external forces act on the system. Mo- 
mentum is a vector quantity so that we may apply the principle only in one 
direction at a time. 

Unlike momentum, we can rarely assume that the energy of a mechanical 
system is conserved. Energy, a scalar quantity, occurs in many forms; in 
mechanics, the common forms are kinetic, potential and strain. In any 
analysis of a system based on energy, possible energy conversion into heat and 
other forms must be taken into account. Such conversion was seen to take 
place in Experiment 7. 


Theory 


For a system of finite bodies acted upon by a single force F, we may write: 
F = X (mb) (8.1) 


where m is the mass of one of the bodies and v is the velocity of the same 
body measured in the same direction as F. Equation (8.1) can be rewritten 
as follows: 


fF dt = Sfm dv (8.2) 


when it is known as the impulse-momentum equation. The time integral of 
the force is the impulse, and this is seen to equal the change of momentum 
of the system. When Fis zero, this leads at once to the principle of the 
conservation of momentum. 

Equation (8.1) can be re-written alternatively as 


fF ds = XJfmv dv (8.3) 


where s is the distance through which the force acts in the direction of the 
force. In this form, we have the work-energy equation. The left-hand side is 
the work done on the system by the external force, and the right-hand side 
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is the increase in kinetic energy of the system. When F is zero, the kinetic 
energy is conserved. 

These equations have analogous forms when applied to rotational motion 
about a fixed axis. We simply replace F, m, s and v as in the table below: 


Linear motion Rotational motion 
force F torque T 

mass m moment of inertia J 
displacement s angular displacement 0 
velocity v angular velocity w 


momentum mv angular momentum Jw 


It is important to measure the relevant quantities about the axis of rotation 
of the system. We can construe the angular momentum to mean the aggregate 
of the moments of linear momentum of the system of bodies, when the term 
moment of momentum is appropriate. 


Experimental Details 


Figure 8.1 shows the important items of a suitable test apparatus. The 
turntable is about 30 cm in diameter and the periphery is flanged to hold a 


Suspension cord 






Cord fixing 
device: 


Ficure 8.1 General arrangement of test equipment. 


strip of recording paper 2 cm wide. It is pivoted as freely as possible by using a 
single steel ball to support the plain end of the shaft, the latter being of 
polished steel of diameter 1 cm; lateral support is provided by a nylon 
bush. A suitable weight for the turntable is 50 N. 

A flat disc, of the same diameter as the turntable and weighing about 
15 N, is supported by a light strong cord attached to the eye of a thin shaft 
fixed to the disc centre. By supporting the disc at a point well above its mass 
centre, it is easy to ensure that when the disc is held stationary its surface 
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is parallel to the surface of the turntable. The rigid supporting frame is fixed 
to the base of the apparatus so that the disc and turntable shafts are coaxial. 
The cord is fastened off by means of a spring-loaded split washer such as is 
used to tension the upper thread on a sewing machine. 

A spigot on the disc is used to locate additional plain discs in order to 
increase the inertia when required; a suitable number of additional discs 
is three, each weighing about 15 N. 

The turntable is put into motion by transferring to it the energy of a 
compressed coil spring. A simple trigger mechanism is used to start the mo- 
tion. The proportions of the mechanism are such that the initial speed 
imparted to the turntable is between 2 and 3 rad/s. A stop is fitted to halt the 
motion after rather less than one revolution. 

Measurement of velocity is carried out by using a freely-vibrating canti- 
lever to make a sinusoidal trace on recording paper wrapped around the 


Ficure 8.2 Measurements from the recording. 


periphery of the turntable. A suitable natural frequency for the cantilever is 
8 to 10 Hz (c/s). It is essential that the pen causes a minimum of retardation 
of the turntable. For this reason a very light stylus is used in conjunction 
with electrically-sensitive paper; an a.c. mains supply is suitable with a 10 kQ 
5 W resistance in series. Pen and bearing friction reduce the initial turntable 
speed by no more than 10 per cent over the whole of the working displace- 
ment (a bigger reduction than this would render the apparatus unsuitable). 
The stylus is released by the same triggering mechanism which initiates the 
turntable motion. 

In use, the horizontal portion of the disc suspension cord is severed when 
the turntable has rotated through between one quarter and one half of a 
revolution. Slipping occurs over about 5° of rotation, after which the disc 
and turntable move together. The velocities of the turntable before the disc is 
released and after slipping has ceased are assessed from measurements of 
the sinusoidal trace on the recording paper. 

Measure the natural frequency of the cantilever by means of an accurately- 
calibrated stroboscope and mark the cantilever accordingly. Measure the 
moment of inertia of the turntable and of each disc by suspending each as a 
trifilar suspension (as described in Experiment 3) and mark the items with the 
appropriate values. 

Carry out the test proper by starting the unloaded turntable and releasing 
the main disc after the turntable has moved through one quarter to one half 
of a revolution. The initial clearance between the disc and the turntable 
is not critical and need be no more than 1 cm. You will have recorded a 
sinusoidal trace as shown in Figure 8.2. With the stylus in its equilibrium 
position, mark the centre line on this recording, and measure AB, BC, CD, 


~ 
N 


MOMENTUM AND ENERGY EXPERIMENT 43 


etc. Plot these measurements against time, as in Figure 8.3, taking the time 
interval between points to be equal to the known period of the cantilever. 
Your plot will represent angular velocity to an appropriate scale, and there 
will be a noticeable step over the interval during which impact occurred. 
Draw the best straight lines through several points in each group either side 
of the step. 

The above lines will not be horizontal because friction will have decelerated 
the system. Neither are they likely to be parallel because the magnitude of the 
deceleration probably will have changed. It is necessary to make accurate 





Angular displacement and velocity 





© Angular displacement 






© Angular velocity 


Ficure 8.3 Angular displacement and velocity. 


assessments of the velocity just before impact and just after slipping ceased 
and, for this purpose, an estimate is required of the time at which impact may 
be supposed to have occurred. Hence now measure AB, AC, AD, etc. from 
your trace and plot these values against time to represent the angular dis- 
placement of the turntable, as in Figure 8.3. Again you will have points in 
two distinct groups. Draw the best straight lines through several points in 
each group adjacent to the interval over which impact occurred. Suppose that 
the point of intersection of this pair of lines represents the time of impact. 
Measure the velocities w, and 2 from the velocity lines corresponding to 
the time of impact as found above. You may suppose that these values 
represent the velocity of the turntable before and after impact. 

You are now in a position to check that, over the short period covering 
the impact between the two bodies, the principle of momentum has been 
obeyed. You may also calculate the loss of kinetic energy due to impact as a 
percentage of the kinetic energy of the turntable just prior to impact. 

Repeat with the main disc loaded with one, two and three additional discs. 
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Discussion 


Did your results confirm that the momentum principle was obeyed to the 
same degree of accuracy in all four tests? If not, discuss the main sources of 
error under the following headings, and suggest possible remedies: 


(a) design of the apparatus, 
(b) execution of the test, and 
(c) interpretation of the results. 


Are you satisfied that the amount of energy lost varied in the expected 
manner according to the amount of the added inertia? To what extent does the 
initial clearance between the disc and the turntable affect the application of 
the energy principle? 
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Forced Vibration of a 
Simple System 


We consider in this experiment the behaviour of a system with one degree of 
freedom as a result of imposing a continuously-varying disturbance from 
outside the system. To simplify the analysis, it is usual to take a disturbance 
which varies sinusoidally. Since the frequency of the disturbance can be 
varied over a wide range, this procedure can provide much useful information 
about the system. 

The system chosen for experiment has the same mathematical form 
as that used in Experiment 6, but the components are arranged in a different 
manner in order to facilitate the tests. 


Theory 


Consider the rotary system shown in Figure 9.1, in which the elastic member 
is a torsion spring of angular stiffness k, one end of which is displaced for- 
cibly through an angle 6, by some external agency; we can refer to this as 






Cc 
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Ficure 9.1 System with one degree of freedom. 


the input to the system. The other end of the spring is attached to and rotates 
with a massive disc which moves through an angle 62, the system output. 
Movement of the output shaft is resisted by a viscous damper of coefficient 
c and the total inertia associated with the output shaft is /. Instead of the 
homogeneous equation (6.3) we now have 


Ib, + cb, + k0, =. ko, (9.1) 
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It is evident that the complementary function solution in 0. of equation 
(9.1) is identical to the solution of (6.3) in 6. In other words, when the input 
shaft is held stationary, the free motion of the output shaft is of the form 
given by equations (6.4) to (6.9), depending on the relative values of k, 
I and c. 

For the purpose of this experiment we are also interested in the particular 
integral solution of (9.1) when 0, = 6, sin wt, i.e. when the input shaft is 
given simple harmonic motion of amplitude 6, and frequency w/2z. 

Because the system is of second order, it must be possible to describe 
it by no more than two parameters. To simplify the analysis it is therefore 
convenient to replace k, J and c by 


and re poe gs (9.2) 


where 7, which has the dimensions of time, is called the time constant and Gs 
which is dimensionless, is called the damping ratio. The latter name arises 
since ¢ is the ratio of the actual value of c to the value it would need to 
assume in order to impose critical damping. Note that 7 is inversely pro- 
portional to the undamped natural frequency. A further useful substitution 
is to put 


A= or = o/b, 


where A is clearly dimensionless and is called the Srequency ratio, i.e. the ratio 
of the forcing frequency to the undamped natural frequency. 

The particular integral solution to equation (9.1) for the given input 
function may now be expressed as 


6. = 0, sin (wt + p)/{(L — A?)? + 42242}4 (9.3) 
where 
= —arctan {2CA/(1 — A?)} (9.4) 


The complete solution in 0, includes, of course, both complementary 
function and particular integral components. In all practical cases the former 
is transitory and soon disappears. In forced vibration work we are concerned 
solely with the latter, called the steady-state solution. Equation (9.3) shows 
this to be an oscillation of the same frequency as the input, having an ampli- 
tude depending on 6,, ¢ and A, i.e. on 6, k, I, c and o. It lags behind the input 
because the phase angle y, which depends on k, J, c and @, is negative. 

Note that, because the system is linear, the output amplitude is directly 
proportional to the input amplitude 6,; we can therefore make 0, any angle 
convenient to the experiment. It is best to express the steady-state forced 
behaviour by the so-called magnification ratio M, meaning the ratio of the 
output magnitude to the input magnitude, together with the phase difference 
y. Hence (9.3) can be replaced by 


M = 1/{(1 — 42)? + 4¢242}4 (9.5) 


FORCED VIBRATION OF A SIMPLE SYSTEM 47 


In Figure 9.2 curves of M and y against 4 for various ¢ are drawn. The 
curves indicate that M for certain ¢ passes through a maximum at a particular 
A. This maximum can be found by differentiating the r.h.s. of (9.5) with 
respect to A and equating to zero. Thus 


Mmax = 1/{22(01 — 22) at 4=(1 — 22%) (9.6) 





Ficure 9.2 Amplitude ratio and phase spectrum. 


When this maximum occurs at 21> 0, we describe the phenomenon as 
resonance. Further analysis shows that resonance can be avoided by making 
£ > 0-707 and equation (9.6) applies only when ¢ < 0-707. Resonance can be 
very destructive in certain engineering applications and it is important that 
the conditions under which it might occur can be foreseen. Note that M — 0 
as A— oo for all €. 


Experimental Details 


It is desirable that the input shaft is displaced with undistorted simple har- 
monic motion. To avoid cyclic variation of o we may employ a Velodyne, 
though, if a relatively large motor is used (0:5 kW), ordinary armature con- 
trol of a d.c. machine would be adequate. The motor drives a crank-slider 
mechanism at about one-tenth motor speed; if space prevents the use of a 
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large connecting rod to crank ratio, a Scotch-yoke mechanism can be 
employed instead. 

The output is converted to oscillatory motion of the input shaft by means 
of a rack engaging a gearwheel keyed to the shaft; a suitable value for 0, 
is 15°. The mechanism should be free from backlash throughout. A suitable 
range of the input frequency /27 is from zero to 2-5 Hz; a speed indicator 
on the motor shaft enables w to be measured. Choosing a spring of stiffness 
about 0-08 N m/rad and a total inertia at the output shaft of about 0-0012 
kg m? gives us a value for the undamped natural frequency of about 1-3 Hz, 
which is mid-way within the range of «/27. 

In this application it is convenient to use a damper in the form of an eddy- 
current brake. The inertia disc may be supported on the shaft of a small 
induction motor which is supplied with direct current; this provides a braking 
torque very nearly proportional to angular speed, as is required for viscous 
damping. In a typical installation the current can be varied to alter the 
effective c, the maximum of 1 A providing rather more than critical damping. 

To record 6, and 0, it is only necessary to fit light pulleys to the input 
and output shafts which drive light cords attached to the pens of a strip- 
chart recorder. The steady-state result would then appear as a pair of sinu- 
soids for which M and y could be determined. This determination is laborious, 
however, and it is not possible to find y with acceptable accuracy. Hence the 
alternative arrangement described below is advisable. 

An oscillating platform replaces the strip chart of the recorder; this is 
driven in a direction perpendicular to that of the output pen by a second 
drive taken from the input motor. The input pen is no longer required. By 
using a differential, a deliberate phase lag is introduced between the motion 
of the platform and of 6,; this angle can be varied by means of a handwheel 
and is indicated by a pointer on a dial. In use, this ‘bias’ lag is made equal to 





Ficure 9.3 Lissajous figure on recording paper. 


y in the system and the value of y appears on the indicator. The record on the 
paper attached to the platform now appears as a Lissajous figure which, in 
general, will have an elliptical shape. When the handwheel is turned until the 
ellipse is converted to a line, the system lag yp is indicated on the dial; ycan 
be determined to an accuracy of +1°, which is some five times better than 
when using strip-chart recording. A typical record is illustrated in Figure 9.3; 
this was made with a light electric stylus on electrically-sensitive paper. 

Of course M is directly proportional to the absolute travel of the stylus 
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since ti is kept constant; it is thus an easy matter to determine M from the 
record. 

In order to work in terms of A, it is necessary first to determine b,. With 
zero damping current, turn the disc through about 180° and release it; time 
the nearly-free motion over about ten cycles. Repeat to check consistency. 

Now set a high value of damping current and switch on the input motor. 
Wary over its entire range and observe the motion of 6,; check the calibra- 
tion of the speed indicator by measuring the period of the oscillations. 

Continue the experiment with less damping; now care must be exercised 
in the region of resonance to avoid damaging the apparatus. Assess the 
minimum value of damping current which should be set when the system is 
taken through resonance. 

The object of the experiment is to produce a family of experimental 
curves to compare with Figure 9.2. Aim at investigating the system behaviour 
for about five constant settings of damping current over the range provided. 
For each setting, carry out two or three tests with 4 <1 and a similar 
number with 4 > 1, making use of the full speed-range of the input motor. 
It will be found that the steady-state condition will always be reached very 
quickly. All that is then necessary is to vary the bias angle to the platform 
drive until the best line recording is obtained. Then read the speed indicator; 
finally read y and determine M from the record. 

Plot M and y against A and sketch in the best curves for each value of 
damping current. 


Discussion 


Make sure that you can obtain equation (9.1) and can solve for the steady- 
state output motion given by (9.3) and (9.4). Differentiate equation (9.5) 
and obtain (9.6) for € < 0-707. 

Try to relate damping current to effective ¢ and also to effective c if I 
and k are known. Does your experimental family of curves differ signifi- 
cantly from Figure 9.2? If so, can you trace the reason or reasons? If, for 
various 2, ¢ for the same damping current changes markedly, this will 
suggest non-linear damping. 

Were any of your line recordings other than straight ? If so, can you suggest 
why? 
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Free Vibration of a System with 
Two Degrees of Freedom 


S 


We now study a conservative system having two principal natural modes and 
an infinite number of combined modes. The experiment introduces the 
phenomenon of beating which occurs commonly in engineering assemblies 
subjected simultaneously to two vibrations of almost identical frequency. 


Theory 


Consider the system shown in Figure 10.1, which consists of two simple 
pendulums coupled together with a spring. For simplicity the pendulums are 





FicurE 10.1 System comprising coupled pendulums. 


similar, each having a length / and a bob of mass m. Of course, they both 
swing in the same vertical plane from supports which are at the same height. 
The spring, of stiffness k, is attached at /, from each point of support and is 
unstressed when the angles of swing 0, and 6, are both zero. There is no 
damper. 

Only the linear behaviour is considered, i.e. the amplitudes are taken 
to be small enough to assume that sin 0, =~ 6, and that sin 6, = 0. The 
equations of motion for the two pendulums may then be written 


ml?6, + mgl0, + k12(6, — 6.) = 0 (10.1) 
ml?6. + mgl0. + k12(6. — 6,) = 0 (10.2) 


TWO-DEGREES-OF-FREEDOM SYSTEM | 


This is a pair of simultaneous second-order homogeneous equations. For 
convenience we introduce the following: 


i,=p+q; %=p-—q (10.3) 


Equations (10.1) and (10.2) may now be written in terms of p and q and, after 
considerable manipulation, we obtain the following independent solutions 
»for 0, and 0.: 


6, = pcos (bit + y1) + cos (bot + po) (10.4) 
8, = pcos (byt + yx) — 4008 (bat + 2) (10.5) 


where p and g are the amplitudes of p and q respectively, b, and y, are 
associated with the motion p and likewise b, and y, with the motion q; 
b, and by are in fact given by 


by = (g/Di (10.6) 
bo = {(ngl + 2KI2)[ml?}4 (10.7) 


In any solution of (10.4) and (10.5), 5, g, yi, and y2 depend on the initial 
conditions. We shall consider several special cases. 


(a) At t= 0, suppose that 6, = 6, = A, 6, = 6, =0. Substitution in 
(10.4) and (10.5) and their derivatives results in y, = yo = 0, p = A 
and g = 0. Hence 


6, = 0, = A cos Dit (10.8) 


Equation (10.8) describes the first principal mode of the system, in 
which the two pendulums have exactly the same motion. The frequency 
is determined by (10.6) which does not contain k. In fact the spring 
remains undeflected during the motion, each pendulum moving in- 
dependently under gravity as if the spring was not present. 


(b) At t=0, suppose that 0, = A, 0.= —A, 6,= 6, =0. Similar 
substitution results in y, = y2 = 0, p = 0 and g = A. Hence 


6, = A cos bet 


9, =~ A cos bat (10.9) 


Now the two pendulums move in anti-phase but otherwise similarly 
at a higher frequency as given by (10.7). This is the second principal 
mode, the spring deflecting so that its mid-point remains stationary. 
For large k or ,, b2 is considerably greater than b, and we say that the 
coupling is strong. For small k or /,, bz is slightly greater than b, and 
we refer to the coupling as being weak. 

We have an infinite number of further choices of initial conditions, 
all of which result in more complex modes than the above, modes 
which combine oscillations of both frequencies b,/27 and b,/2z. 
We choose a third special case. 
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(c) At t=0, suppose 0, = A, 6, = 6, = 62 =0. Similar substitution 
results in yp, = p2 = 0, p = 4A and g = —}A. Hence © 


6, = —4A(cos bet — cos byt) (10.10) 
6. = $A(cos bot + cos bjt) (10.11) 
which may be written as 
6, = {A sin 4(b. — b,)t} sin (bo + 5i)t (10.12) 
6, = {A cos 4(bz — b,)t} cos (bg + b;)t (10.13) 





Ficure 10.2 The phenomenon of beating. 


We can translate equation (10.12) as meaning an oscillation of fre- 
quency mid-way between the principal mode frequencies, having an 
amplitude which varies sinusoidally. For a weak coupling, the varia- 
tion of amplitude is slow. The motion is illustrated in Figure 10.2 
(upper part) and is an example of beating. The period at which the beat 
recurs is 27/(b, — b,) and hence the beat frequency is (6, — b,)/2z. 


Beating becomes evident whenever a structure is simultaneously subjected 
to two vibrations of similar frequency. For the special case above, where the 
two vibrations are of the same amplitude, see equation (10.10), the resulting 
amplitude becomes zero once every beat period. 

We note that, for our present system, the two pendulums vibrate in the 
same manner, see (10.12) and (10.13), except that they beat alternately; 


~~ 
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see Figure 10.2. Being a conservative system, the total energy remains con- 
stant. As the vibration of one pendulum dies away, that of the other builds 
up due to the transfer of energy from one part of the system to the other via 
the spring. 


Experimental Details 


It is desirable first to match a pair of pendulums, about 80 cm long, and to 
ensure that each is suspended in bearings as free from friction as possible. 
A spring should be chosen and fitted such that 2k/? is between 0-25 and 0:35 
of mgl to provide a weak coupling. To ensure that the spring retains the same 
stiffness in compression as in tension, it should be made with a small number 
of turns of large diameter and anchored to the pendula in bearings free from 
friction and from backlash. Of course the spring should be unstressed when 
the system is in equilibrium. 

Remove the spring and set the separate pendulums in motion in phase 
with each other with a similar amplitude not exceeding 10°. Check that 
they remain in phase over a long period; otherwise lengthen one of the 
pendulums accordingly. Also check that both decay in amplitude at about 
the same rate; it is desirable that this rate of decay is as small as possible 
though it cannot be reduced to zero in practice because of friction and 
windage. 

Re-fit the coupling spring and check that it is unstressed when both pendu- 
lums hang freely. Practise setting the system in motion in its principal modes; 
this is not a simple matter because of the special nature of the modes. Measure 
the two periods over a large number of cycles and calculate b, and bp. 

Now set the system to beat. This is comparatively simple since it is necessary 
only to release one pendulum from a displaced position at the same time as the 
other is released from its equilibrium position. Measure the beat frequency 
and check its relationship to b, and b2. Measure the frequency of swing of 
each pendulum and check whether this is given by the mean of b; and by. 
Study carefully the overall motion of the system and check that the pendulums 
beat alternately. Also check whether the amplitude of each becomes zero 
during each beat period; if not, you have not started the motion precisely 
according to the conditions stated in (c) above. Check that the pendulums 
swing 90° out of phase with each other, i.e. that one is at maximum displace- 
ment when the other passes through its equilibrium position, as suggested 
by equations (10.12) and (10.13) and by Figure 10.2. 

Now set the system in motion without attempting to impose any special 
initial conditions. In general you should still find that each pendulum 
beats but that the amplitude of the beat no longer reduces to zero. 


Discussion 


Check that you are able to write down equations (10.1) and (10.2), and that 
you can obtain (10.6), (10.7) and (10.8) to (10.13). 
Describe the ease or otherwise of initiating the various kinds of motion 
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called for in the procedure; what conclusions do you draw from your 
experience ? 

To what extent did frictional effects make the experiment difficult to 
carry out or the various frequency components difficult to measure? Were 
you able to assign causes to the various ways in which the experimental 
behaviour departed from the theory? 

Describe examples of beating which you have experienced, e.g. (a) the 
vibration of a diesel railcar when stationary with the engine idling, (b) the 
motion of the floor of the turbine hall of a power station, (c) piano tuning. 

Analyse the process of energy transfer during beating between the two 
pendulums of the system you have tested; you should start with an expression 
for the energy in terms of A, g and the system parameters. 


11 
Shaft Whirling 


Whenever a shaft to which rotors are attached is caused to rotate, it is found 
that whirling occurs at a certain speed or speeds however well the system is 
balanced. This means that the shaft has deflected and rotates in its bent form. 
Often whirling does not occur within the normal operating speed range of a 
system, but the phenomenon is common whenever a fairly flexible shaft 
and/or several heavy rotors are involved, e.g. turbine machinery. Because of 
the excessive stresses involved, a knowledge of the whirling speeds of a system 
is essential in order that such speeds are avoided when the system is operated. 

This experiment investigates the phenomenon for a system with one 
degree of freedom, and shows how it is related to the transverse vibration 
of the same system. 


Theory 


First consider the free transverse vibration of a simply-supported uniform, 
light, straight beam, of length /, to which a single concentrated body, of mass 
m, is attached at distance p from one end, as in Figure 11.1. 


Mass m 
p 


fe 








Ficure 11.1 Transverse vibration of a beam. 


Elementary bending theory (see Experiment 21) informs us that the force F 
required at the body to deflect the body by 0 is given by 


F = 3EII6/(pq)? (11.1) 


3—11 
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where E and Jare respectively Young’s modulus for the beam material and the 
relevant second moment of area of the beam section. Hence, in transverse 
vibration, the beam may be assumed to act as a simple spring of stiffness k 
given by 


k = Fé = 3Ell|(pq)" (11.2) 


By comparison with equations (6.6) and (6.7), the frequency of free transverse 
vibrations is given by 


bn _ (k[m)! 
Qn = arr 
Ell|m)* 
where b= eee) (11.3) 
PY 


Note that we have ignored damping from any cause and beam deflection 
due to shear; also that, by neglecting the mass of the beam, we obtain only 
one mode of vibration. Clearly, as the body is moved along the beam, the 





p 


FicureE 11.2 Variation of natural frequency with the position of the body. 


natural frequency changes, being a minimum when the body is sited mid-way 
between the ends; this variation of frequency is illustrated in Figure 11.2. 
We have ignored the tilting of the body when sited away from the mid-point. 

Now consider a uniform, light, circular shaft of length /, having the 
same E and J as above, supported freely in self-aligning bearings at its ends, 
as in Figure 11.3. A disc of mass m is firmly attached to the shaft at distance 
p from one bearing, and the whole is rotating at angular speed w. Figure 11.3 
also shows a view in the direction from one bearing to the other, in which 
O is the centre of rotation (the centre-line of the bearings), A is the centre 
of the deflected shaft (the geometrical centre of the disc) and C is the centre of 
mass of the disc. The distance OA is r and the distance AC is the eccentricity e. 
In practice e is often very small but it can never be assumed to be zero; 
in general AC is not in line with OA. The configuration OAC remains fixed 
during rotation and the shaft is said to whirl; although there is no change in 
the deflection of the beam, nevertheless it is convenient to assume that some 
aerodynamic damping force exists and, furthermore, to assume that this 
acts radially inwards, i.e. in the same direction as the force due to the beam 
stiffness in bending. 
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Ficure 11.3 Whirling of a shaft. 


Taking the co-ordinates of A as x, y, we may write the equations of motion 
for the displacements of the body in the X and Y directions as follows: 


mD*(x + e cos wt) + cDx + kx = 0 (11.4) 
mD*(y + e sin wt) + cDy + ky =0 i 


where c is the ‘damping coefficient’ and D is the differential operator d/d¢. 
The solutions of the above in the steady state are 


x =rcos (wt + y) 

y=rsin(ot + y) he) 
where r= el?/{((l — A)? + Alea}? (11.6) 
and y = —arctan [2CA/(1 — 4*)] 


where A = w/b, and ¢ = c/2(mk)!, these symbols having been introduced 
in Experiment 9. 

Ignoring the damping, the shaft deflection at the body, 7, becomes 
e22/(1 — A?). This shows that when w = b,, i.e. when A = 1, r is infinite. 
This is the condition which must be avoided in practice to prevent violent 
deflection and sudden fracture of the shaft. The critical speed is b,/2m 
and is identical to the frequency of free transverse vibrations of the shaft. 
The phase angle p is zero when 2 < 1 and 7 when A> 1; this means that 
O, A and C are in line when damping is absent. Below the critical speed, 
A is between O and C, and above the critical speed C is between O and A. 
For 2> 1, C approaches O more nearly as the speed is increased, which 
explains the smoothness with which a shaft runs above its critical speed. 

The critical speed varies with the position of the body along the shaft 
in the same manner as the transverse mode frequency, as shown in Figure 
wt i2! 

When damping is not ignored, it is seen that p> 0 for A — 0 and that y 
increases in magnitude gradually to 47 as A is increased to unity. For Al, 
the magnitude of y increases further and reaches 7 at A= ©. 
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In the above analysis, we have ignored shear deflection, the mass of 
the shaft and the gyroscopic effect (see Experiment 12) due to the tilt of the 
disc when situated away from the mid-point of the shaft. The latter increases 
or decreases the critical speed depending on whether the gyroscopic precession 
is in the same or opposite direction compared with the direction of rotation 
of the shaft. 

It must be appreciated that when more than one disc is attached to the 
shaft and/or when the shaft mass cannot be ignored, there will be more than 
one critical speed. 


Experimental Details 


A suitable apparatus is shown in Figure 11.4. This utilizes a steel shaft, of 
diameter 1 cm, supported horizontally in two spherical self-aligning bearings 
80 cm apart. A pair of guard rings is fitted to prevent excessive deflection of 
the shaft. A disc, of diameter 10 cm and mass 1-2 kg, is fitted with a collar 
to enable it to be fixed at any point along the shaft. These values result in a 
natural transverse mode frequency, by equation (11.3), of about 14 Hz 
when the disc is situated at the mid-point of the shaft. 


Application i 
of shaker Guard Guard Coupling Shaker 











Ficure 11.4 Test apparatus. 


For whirling tests, one end of the shaft is driven, through a flexible coupling, 
by a variable-speed electric motor having a maximum speed of 2000 rev/min. 
A radial line marked on the disc enables a stroboscope to be used to measure 
the shaft speed. 

For transverse vibration tests, an electromagnetic shaker is used to 
oscillate the shaft horizontally at a point close to one of the bearings. The 
shaker is driven by a power oscillator, the combination being capable of 
providing a maximum force amplitude of 10 N, the displacement amplitude 
being restricted to 2mm, over a frequency range from 5 to 10 000 Hz; the 
upper end of this range is not required in the test. 

With the shaker disconnected, fix the disc mid-way along the shaft and 
start the electric motor. Slowly increase the motor speed until the critical 
speed is reached, i.e. until the deflection reaches a maximum (damping 
provided by the bearings should obviate excessive deflection). Increase 
the speed beyond critical and note the smooth running of the shaft. Return to 
the critical speed and measure this by means of the stroboscope. 
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Fix the disc at a minimum of ten different measured positions along the 
shaft and find the critical speed in each case. 

Connect the shaker, return the disc to mid-way along the shaft and apply 
a small-amplitude force at the lowest possible frequency. Slowly increase 
the frequency until resonance is obtained, i.e. until the natural transverse 
mode is excited. Note the frequency on the calibrated dial of the oscillator. 

* Repeat for the same positions of the disc as used in the whirling tests. 

Calculate b,/27, using equation (11.3), for the range of values of p used 
in the tests. 

Plot the theoretical frequency against distance of the disc along the shaft. 
Also plot the critical frequency values found in the whirling tests and the 
natural frequency values found in the transverse vibration tests to the same 
axes. 


Discussion 


Compare the three curves you have plotted and discuss the discrepancies 
between them, taking into account the following: 


(a) mass of the shaft, 

(b) shear deflection of the shaft, 

(c) damping of the motion, 

(d) twisting of the disc in transverse vibration, and 
(e) gyroscopic effect in whirling. 


12 ) 


Simple Gyroscope 





The precessional behaviour of a gyroscope, as described in an analytical 
treatment, is somewhat surprising when met for the first time. It is therefore 
particularly valuable to examine this behaviour in the laboratory. 


Theory 


Referring to Figure 12.1, consider a thin circular disc attached to a shaft 
lying in the axis Z,Z,. The instantaneous plane of rotation of the disc is 





Y2 


FicureE 12.1 Gyroscopic precession. 


X,0' Y; and the disc rotates at uniform speed «, in the direction indicated. 
At the same time the shaft Z,Z. rotates about some vertical axis Y, Y, at 
uniform speed w, in the direction shown; by this means the disc is said to 
precess. 

Consider a particle P in the disc, of mass dm, rotating about O’ at radius r. 
At some instant, O’P makes an angle 6 with O’X,. The velocity x of P 
relative to O’ in the direction X;X4 is given by 


x = w,r sin 0 
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But Z,Z, is rotating about Y; Y, at speed w.; hence 
Coriolis acceleration of P = 2m, (ar sin 0) = 2m,w.y 
This is in the direction Z,Z,. It follows that 
Inertia force on P = 2m,.dmy 


which is also in the direction Z.Z,. This force has a clockwise torque about 
X2X, given by 


Inertia torque = 2,w,dmy? 


For the entire disc, we may write 


Total inertia torque = 2w,w.Ldmy? 
= 20, Wo, 


where J, is the moment of inertia of the disc about its diameter X3X,. But 
I, = 4, where J is the polar moment of inertia about Z,Z.. Hence 


Total inertia torque = w,@2J (12.1) 


Suppose that the shaft Z,Z, is supported freely at O; then the above 
torque can be provided by hanging a body of mass M at Q, distant z from 
O, such that 


Megz => G = O1WoT (12.2) 


The above analysis shows that when a couple C is applied to the system, 
it does not turn about the axis of the couple but about an axis which is 
perpendicular to both the couple axis and the axis of spin of the disc. Of 
course, there is a rule connecting the directions of 1, 2 and C; a reversal 
of couple direction requires the reversal either of w, or of 2 and the couple 
can remain unchanged if w, and w. both change. Note that Y; Y. can be any 
vertical axis including Y3 Y4. 

For a full examination of gyroscopic precession, the reader is referred to 
Arnold and Maunder [1]. 


Experimental Details 


Suitable apparatus is shown in Figure 12.2. A straight bar supports an 
electric motor driving a circular disc at one end and a balance body at the 
other end. The speed of the motor can be varied in the range 1000 to 6000 
rev/min. The bar is supported at its centre of mass on a rigid stand by means 
ofa cup-and-cone pivot bearing. Current is fed to the motor through a mercury 
cup surrounding the bearing into which dips a stainless steel pin attached to 
the bar; the current is returned via the pivot bearing. 

Grooves are cut in the bar at 2-5, 5-0, 7:5, etc. cm from the pivot to locate a 
hanger for supporting one or more small bodies. The gravitational force on 
the body, together with the pivot reaction, provides the applied couple C. 
A stroboscope is used to illuminate the disc in order to measure ,. The 
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bar precesses slowly about its pivot, and to measure wz, it is only necessary 
to provide means for assessing accurately the time for a few complete 
revolutions. 

Set the motor speed control at maximum and allow several minutes for 
the speed «, to settle. Hang a small body in a groove near to the pivot and 





Ficure 12.2 Simple gyroscope apparatus. 


observe the slow precession of the bar. Now hang the body on the other side 
of the pivot and note that the precession reverses its direction. Should the 
bar gradually assume an angle of inclination to the horizontal, it must be 
returned to the horizontal by applying light finger pressure horizontally 
on the appropriate side of the pivot. Care should be take not to set the 
bar into oscillation about the pivot. 

For the test proper, run the disc at about four different values of @. 
In each case, apply various couples in both directions to the bar by using the 
small bodies of known mass. In a typical apparatus having a disc of inertia 
about 2 x 10-* kg m?, a suitable maximum value for the applied couple C 
is 0-7 Nm. Measure w. and check @. 

Plot w, against C for each w,. If J is known, draw the theoretical curves 
on the same diagram by using equation (12.2). 


Discussion 
Compare the test points with the theoretical curves and try to explain any 
discrepancies. 

What is the effect of the small frictional drag at the pivot? Will this have 
affected the accuracy of the experiment? 
Reference 


1 .ARNOLD,R.N.,and MAUNDER, L., Gyrodynamics and Its Engineering Applications, 
Academic Press, 1961, 
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First-Order Process Lag 





This very common dynamic process lag is met in many branches of engineer- 
ing [1]. A knowledge of its behaviour is essential in the application of linear 
control system theory. 


Theory 


A process having an input 6, and an output 6,, both functions of time and 
preferably having the same dimensions, may obey the following linear 
equation: 


76, + 0, = 0; (13.1) 





FiGureE 13.1 Liquid level process. 


where 7 has the dimensions of time and is called the time constant. Equation 
(13.1) may be expressed more conveniently in terms of the differential 
operator D; thus 


1 
6.= (<7) 6, = T(D)6, (13.2) 
where the so-called transfer operator T(D) relates the output to the input. The 
transfer operator in this case is 1/(1 + 7D) and is known as the exponential 
lag or the first-order lag. 

The process chosen for experiment is shown in Figure 13.1. Incompressible 
liquid is contained in a pair of tanks connected by a restriction. The gravita- 
tional head 6, in the left-hand tank is regarded as the input, and likewise 
6. is the output. It is assumed that the flow through the restriction is laminar. 
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Hence, if J is the volume flow-rate and R is the resistance to flow, we may 
write 

6, — 0, = RI (13.3) 
Thus R is a property of the restriction, meaning the volume flow-rate/pressure 
drop, and depends on (a) the physical dimensions of the restriction and (b) 


the properties of the liquid. The flow through the restriction must equal the 
flow into the right-hand tank; hence 


I= CD60, (13.4) 


where the capacitance of the process C is in this case simply the cross-sectional 
area of the right-hand tank. 





FicureE 13.2 Step response of first-order lag. 


Eliminating J between (13.3) and (13.4) and putting 7 for CR, we obtain 
(13.2). Hence the process is an example of a first-order lag. 

It should be noticed that the cross-sectional area of the left-hand tank 
does not enter the analysis. However, it is essential for C to be constant 
throughout the right-hand tank. ; 

The apparatus available for test lends itself to the measurement of the 
transient behaviour of the process, in particular the transient response in — 
6, to a step change in 6,. It is convenient to measure 6, and 0, from some 
datum level above the restriction to permit of negative values occurring. 
Suppose that initially both are zero. Now suddenly make 0; = 9, and keep 
it constant at this value. The relevant solution of (13.1) or (13.2) may be 
written [1] 


0, = 6,(1 — e-"") (13.5) 


and is illustrated in Figure 13.2 together with the step in 6,. 

Eventually 6. becomes equal to 6, but in general it lags behind 6,; it 
approaches 6, exponentially. This explains the name ‘exponential lag’. 
Two properties of this response curve are worth noting. The initial slope is 
6,/7 and 0, reaches 0-865 6, in time 2r. 
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In practice it is unlikely that the flow through a restriction will be laminar 
except at very low rate. It follows that we cannot expect the linear analysis 
above to be obeyed precisely in the test apparatus where a commercial 
valve is employed. The value of the experiment lies in noting the differences in 
behaviour of a real process compared with the behaviour predicted by the 
simple theory. 


Experimental Details 


Two sections of Perspex cylindrical tube of 8cm diameter are mounted 
with their axes vertical as shown in Figure 13.3. Cylinder A is closed at both 


Recording drum 


Straight-line __ 
linkage 





Float 


Input head scale 


Valve C 





Air 


cece Set al | 
Ficure 13.3 Apparatus for process-lag response tests. 


ends but cylinder B is open at the top. The cylinders are connected by a pipe 
of diameter 1-5 cm containing a commercial screw-down valve C. A mano- 
meter with a scale is fitted to the pipe on the left-hand side of C. The assembly 
is filled with water to the level indicated. 

A hollow brass float in B is connected to a straight-line linkage to move 
a pen vertically over a recording chart fitted to a drum which is driven at 
a uniform speed about a vertical axis by a small electric motor. 

The effective liquid head to the left of C is adjusted by admitting air under 
pressure into A. The air is supplied from a reservoir which is pressurized 
by a hand pump to about 15x 10° N/m’. Air is admitted through the 
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needle regulating valve P and exhausted to atmosphere through a similar 
valve Q. 

With all the valves open, check that the water level is at an appropriate 
point on the manometer scale. Check the recording drum drive and measure 
the paper speed, which should be between 20 and 25 cm/min. Close valve P 
and pressurize the air reservoir. 

Set the screw-down valve C at half-turn open. Close valve Q. Start the 
drum and, when the pen has traversed about 2 cm, gently open valve P and 
regulate it continuously so as to maintain the level in the manometer at 10 cm 
above the starting level. This will cause the float to rise in cylinder B. As 
the float motion terminates it will again be necessary to close valve P 
completely. Stop the drum. 

With P closed, start the drum and gently open valve Q and regulate it so 
as to maintain the manometer level at its original reading. At the termination 
of this second response test, valve Q will be open. Again stop the drum. 

You have tested the process comprising B and C with forward and reverse 
steps of input head. You have recorded curves of output head change against 
time. Carry out similar tests for step changes of 2, 4, 6, 8, and 10 cm, with 
C one-quarter, one-half, three-quarters and one turn open. 

For each part of a selection of recordings, determine 7 by measuring (a) 
the initial slope and (b) the time to reach 0-865 of the step. Under certain 
circumstances these will not be the same because of non-linear behaviour. 
Construct a true exponential curve through the origin based on the 7-value 
found in (b). 

Now take the record for the 10 cm step when the valve was one turn open. 
Measure (6, — 6.) and D6, for a selection of points on each curve. Plot 
logarithms of these values. Note that, in practice, (0; — @2) varies with 
(D6.)" where n is a number other than unity. From your graph note how n 
varies with (0, — 6,). 


Discussion 


How does the real process lag differ in behaviour from a linear exponential 
lag? What changes occur in the apparent ‘time constant’ of the lag as the 
magnitude of the step is altered ? How does the lag vary with valve opening? © 
How does the direction of motion through the valve affect the lag? 

What restriction must be placed upon input changes in order that the 
transfer operator given by equation (13.2) may represent the process? 


Reference 
1. WEBB, C. R., Automatic Control, McGraw-Hill, 1964. 
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Position Control System 


Following the last test which investigated the behaviour of a simple transfer 
lag, this experiment deals with a complete closed-loop control system, 
arranged to control the angular position of a shaft. We start with a basic 
proportional system and then add conventional compensation to improve 
the performance. 


Theory 


Let 6, be the angular position of a shaft to which a disc of inertia / is attached, 
as in Figure 14.1. The motion of the shaft is resisted by viscous friction of 





Torque 6, 


Ficure 14.1 Physical system to be controlled. 


coefficient c and by an external load torque L. In the basic system, torque is 
applied to the shaft by an amplifier and servo motor; this torque is propor- 





Ficure 14.2 Basic block diagram. 


tional to the error 6 between the desired position 6, and the actual positi n 
6,. The block diagram is shown in Figure 14.2, where the transfer operator 
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represents the controller (amplifier) and servo motor and G, the process 
parameters J and c. The following equations apply [1]: 


kO —L = (ID? + cD), (14.1) 
= 0; P= 0, s (14.2) 


where k is the proportional control coefficient, i.e. (input torque)/(position 
error). 
Eliminating 6, from (14.1) and (14.2), we obtain 


UD? + cD + k)6 = UD? + cD)6,; + L (14.3) 


which gives the error as a function of 0; and L. The coefficient of 9 shows the 
system to be mathematically similar to the free and forced simple vibrating 
systems considered earlier. The quantities defined by equation (9.2) apply 
here also; we restrict our attention to the underdamped case, ie. € < 1. 
Whenever a disturbance is made in 6; or in L, the system response is damped 
oscillatory as given by equation (6.5). 

We are also concerned with the error which persists after the transient 
has faded, i.e. the steady-state error 6,,. By equation (14.3), for constant L and 
6,, the steady-state error is given by 


O55 = L|k (14.4) 


Hence, in the absence of external load, 0,, is zero. Note how the controller 
(amplifier) gain, which contributes to k, affects the value of 8,,. 

When the desired speed is suddenly changed from zero to D0; and held 
uniform, i.e. a ramp input function is applied, the steady-state error, now 
called the velocity lag, is given by 


955 = (c/k)DO, + Lik (14.5) 


Clearly, steady-state errors can be reduced by increasing k, and the 
consequent increase in the response frequency b/27 is beneficial. However, ¢€ 
becomes smaller which is unstabilizing, and hence the simple system cannot be 
stable without exhibiting velocity lag. An equally important drawback of the 
simple system is the inherent load imposed on the servo motor by the viscous | 
drag but, of course, the latter is required for stability. 

We first seek a means. of artificially replacing the viscous drag which 
will permit the use of a less powerful controller and servo motor but without 
loss of stability. We measure the speed D9, of the output shaft and apply an 
additional torque —k,D0, to the process; this is called output-velocity 
feedback and results in (14.5) being modified to 


O55 = [(¢ + ka)/kIDO; + L/k (14.6) 
The parameter c is no longer required for stability since k, takes over this 


function. We may therefore specify bearings exerting negligible resistance and 
consequently a less powerful servo motor. 
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We now add a further path to the system to reduce the velocity lag. 
Another additional torque, +k,D6;, called input-velocity feedforward 
(because of its positive sign), is applied to the process. Omitting c altogether, 
the steady-state error now reduces to 


6., = Lk (14.7) 


Thus in the absence of external load, there is no velocity lag. 


Experimental Details 


Small components as used in instrument servo mechanisms are employed 
in the apparatus described below and shown in Figure 14.3, which is one of 
several suitable arrangements. 

A two-phase 115 V 400 Hz a.c. servo motor A drives the output shaft 
through reduction gearing, a scale on the shaft indicating 0,. Viscous friction 
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FicureE 14.3 Layout of components of position control system. 


is provided artificially by an eddy-current brake, the current through the coil, 
and therefore c, being adjustable. A disc is keyed to the motor shaft to 
provide the inertia J. For experimental purposes, I is made large deliberately 
to ensure a low operating frequency in the range 0-3 to 0-5 Hz. The actual 
position 05 is measured by the synchro transmitter B and DO, by the tacho- 
generator C. A band brake acts as a dynamometer and provides the external 
load L. 

The desired-value shaft can be turned by hand or driven at various values 
of uniform speed in either direction by a d.c. motor C. The actual position 
0,, which is indicated on a scale, is measured by the synchro transmitter EF, 
and D6, is measured by the tacho-generator F. For convenience, a mechanical 
differential enables the error 0 to be indicated on a third scale. 
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The synchros supply an a.c. amplifier which in turn supplies the control 
phase of the servo-motor, & being adjustable by a knob on the amplifier. 
Signals from the tacho-generators also feed the amplifier, there being a 
separate knob to vary ky. 

The steps given below are suggested in order to demonstrate the successive 
improvements in the system indicated by the theoretical outline. 


(a) 


(b) 


(c) 


(d) 


(e) 


(f) 


(g) 


(h) 


(i) 


(j) 


First set up the equipment as follows: 

(i) set k, at zero, 

(ii) set k at maximum, 
(iii) set ZL at zero, 
(iv) declutch the drive from D, and 

(v) switch off the eddy-current brake. 
Switch on the amplifier. The equipment should operate according to 
equation (14.3) with c and L both equal to zero, i.e. with 6; constant 
the 6, and @ shafts should oscillate continuously (hunt) because of the 
lack of damping. Measure the undamped natural frequency 5,/27. 
Switch on the eddy-current brake and set c at maximum. Apply a step 
in 6; and note the highly-damped response in 0, and 6; check that 0 
returns to zero. Note also that the maximum speed of the servo motor 
is limited. Try several steps with smaller c. : 
Switch off the brake. Switch on the output-velocity feedback and set 
k, to maximum. Try a step in 6; and note that there is a similar response 
to the above except that the servo motor maximum speed is much 
greater. Repeat with smaller k,. 
Switch off the output-velocity feedback. Apply a moderate load L at 
the dynamometer. Try a step in 0; and verify that Z stabilizes the 
response. Note the steady-state error. For a large L, verify that the 
maximum speed of the motor is limited. 

Re-set L to zero. Connect the drive from the motor D to the 6; shaft. 
Switch on D and adjust D0; to about | rev/min. Note that 0, follows 
in a non-uniform manner and that @ hunts about zero. Reverse the 
direction of D and check that the response is similar. 
Switch on the eddy-current brake and adjust for a moderate value of c. 
Check that the hunting is damped out but that there is a steady 
velocity lag. Also observe that the servo motor can only just overcome 
the viscous drag at the desired speed demanded. 
Switch off the eddy-current brake. Apply the output velocity feedback 
with moderate k,. Observe that the rapid and stable response to the 
ramp demand even when D6; is increased to 2 rev/min and that the 
velocity lag remains. 
Apply input-velocity feedforward. Check that the velocity lag dis- 
appears. Reverse the direction of D and further check that the response 
is rapid and stable and without steady-state error. 
Apply a small load L at the dynamometer. Check that this causes a 
velocity lag which increases with L until the servo motor is no longer 
able to keep up with the demanded speed. 
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Discussion 


Make sure that you can formulate equations (14.1) and (14.2) and that you 


can obtain the derived equations. 
Compare the results of the various tests you have carried out with the 


linear theory given and try to explain any anomalies. 


7 
Reference 
1. WEBB, C. R., Automatic Control, McGraw-Hill, 1964. 
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Transverse Vibration of an 
Elastic Beam 





This experiment is concerned with the vibration of a system with an infinite 
number of degrees of freedom. This occurs whenever the mass is distributed 
instead of being concentrated at one or more points in the system. We shall 
confine our attention to the transverse vibration of a cantilever due solely to 
the mass of the cantilever itself. 


Theory 


Consider a uniform slender straight beam, of length / and mass m, supported 
in some manner which we shall discuss later. Suppose that the beam is 
vibrating transversely and freely. At a distance x along the beam from some 
arbitrary point, let the deflection at some instant be y. Consider a section 
6x long having mass (m/l)dx; this requires an accelerating force F given by 


m 
Fam ox od (15.1) 


This force arises because of the change of shear force (0Q/éx)éx, where 
the shear force for simple bending deflection is given by 


BORK | 
Q=-E=zx 3 (15.2) 
Therefore 
P= 22 by = BF ox (15.3) 
Eliminating F from equations (15.1) and (15.3), we obtain 
arty m Oy 
dxt = TET Oe (15.4) 


We may assume that the motion at x is simple harmonic, i.e. 


y =F cos byt (15.5) 
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where j is the amplitude at x and b,/27 is the natural frequency, damping 
being assumed negligible. Substituting this expression for y in equation 
(15.4) and simplifying we get 


OF 7 Mee 

axt  IEI 
which, for convenience, is re-written 

atp 

oa ee =O (15.6) 
where 

= (a/)?(El/ml?)! (15.7) 


The solution of equation (15.6) is 
¥=Acosax + Bsin ax + Ccosh ax + D sinh ax (15.8) 


In the above equation, A, B, C and D are arbitrary constants depending on 
the end conditions of the beam. Variations are numerous: the beam can be 
freely supported at both ends, fixed at one end and freely supported at the 
other, etc. For the particular case of a uniform cantilever considered here, 
suppose that the origin is at the fixed end. Then, since the deflection and 
the slope are both zero at the fixed end we have 


First mode 





Second mode 





<i Third mode 





———— 0:906 { ———— 
0-644 ber 


Ee ee et 


Ficure 15.1 First four natural modes. 






Fourth mode 
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Also, since the bending moment and the shear force are zero at the free end, 
we have: 
ey ey 
axa = ae 2 atx = 1 
Substituting these four conditions into equation (15.8) and its derivatives 
leads to the condition 


cos al cosh a/ = —1 (15.9) 


Values of a/ satisfying equation (15.9) are 0°87, 4-69, 7°85, 11-0, .. ., 
enabling 5, for the first four modes to be found from equation (15.7). The 
cantilever shapes when vibrating in these modes are shown in Figure 15.1, 
the theoretical positions of the nodes being indicated. 


Experimental Details 


The material and dimensions of the cantilever to be tested should be chosen 
so that the first four natural frequencies are within the range of a single 
exciter; also the cantilever should be sufficiently flexible to enable the dis- 
placement amplitude at the modes to be large enough to enable the mode 
shapes to be clearly discernable with the naked eye. A suitable test cantilever 
is one made of mild steel strip, of free length 7 = 1 m, width b = 12-7 mm 
and thickness d = 3-17 mm. Such a cantilever has a density of 7700 kg/m, so 
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Ficure 15.2 Arrangement of test cantilever. 


that m = 310 g. I about the appropriate neutral axis (parallel to the width 
b) is found to be 33-4 mm‘. Taking E = 207 x 10° N/m?, the theoretical 
value of b, for the first mode may be found from equation (15.7) to be 16-6 
rad/s. This corresponds to a frequency of 2:64 Hz. Similarly frequency values 
for the next three modes may be calculated as 16-5, 46-3 and 90.7 Hz. 

It is convenient to support the cantilever horizontally, and the clamp 
must anchor the fixed end as rigidly as possible. This may be arranged 
so that the cantilever is free to move either in a vertical plane or in a hori- 
zontal plane; the former is preferable to facilitate visual observation, and the 
dimensions suggested are such that the cantilever does not deflect unduly 
under its own weight. Figure 15.2 illustrates the general arrangement. 

A low-power source of excitation of variable frequency may take the 
form of an electromagnetic shaker or a simple mechanical shaker driven by an 
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eccentric. It is important that the excitation is applied at a point which 
interferes least with the free transverse motion of the cantilever, and this 
dictates that the distance a in Figure 15.2 is kept small compared with /. For 
the cantilever suggested, a suitable value for d is 5 cm and, when a positive- 
displacement mechanical shaker is employed, an excitation amplitude of 
about 0:10 mm is suitable. The eccentric may be driven by a high-speed 
electric motor of about 50 W with means for varying the speed sensitively 
from zero to 10 000 rev/min. Safety considerations demand that a guard is 
fitted to restrict the amplitude at the free end of the cantilever to about 2 cm. 

A more elaborate means of excitation may be adopted to reduce the 
discrepancies between the test and theoretical mode frequencies. This uses 
a light inertia-type shaker fixed to the beam and driven by a flexible shaft. 

To determine frequency, it is best to measure the speed of the exciting 
motor with, for example, a magnetic impulse counter. 

The amplitudes at resonance should be sufficient to enable the nodes to be 
detected easily, but the assistance of a stethoscope may be found useful. 

The first test suggested is to set the cantilever in free motion by hand (with 
the shaker removed). Is this the first mode? Check the frequency by a 
sufficiently accurate method. 

Re-connect the shaker and slowly increase the motor speed from zero 
until resonance occurs at the first mode. Inherent damping should prevent the 
amplitude becoming excessive. Check the frequency by measuring the motor 
speed and compare with the value found above. Further increase the motor 
speed until the second mode is excited. Again measure the frequency. Estimate 
the position of the single node by marking the cantilever where the amplitude 
appears to be a minimum and by subsequent measurement when the canti- 
lever is stationary. 

Now excite the third and fourth modes and repeat the above measurements. 
Even for the fourth mode, the amplitude should be sufficient to detect the 
cantilever shape easily by eye. 


Discussion 


Compare the test and theoretical values of natural frequency and the actual 
node positions with the theoretical positions given in Figure 15.1. Suggest 
reasons for the differences found. 
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Vibration of a Multi-mass System 


We now consider systems in which the body masses are concentrated, 
where the bodies are connected by elastic members assumed to be massless 
and where frictional damping is neglected. As in the previous experiment, we 
have a large number of degrees of freedom but now the number is finite. 
Analytical methods are available for such systems which, in many cases, 
reliably predict real behaviour; however, the analysis becomes increasingly 
tedious the greater the number of degrees of freedom. An experimental 
method is useful, therefore, especially when the general pattern of behaviour 
is known from theoretical considerations. ; 


Theory 
Consider the torsional system shown in Figure 16.1, which consists of a 


light elastic shaft to which are attached five rotors. For simplicity, rotors 1 
to 4 each have a moment of inertia of J; rotor 5 has a moment of inertia of 


A sin ot} 


Ficure 16,1 Multi-mass torsional system. 





20 I. Also, the stiffnesses of the shaft sections connecting rotors 1 to 4 are 
each k, while that of the section between rotors 4 and 5 is 0-5 k. There is a 
light extension shaft of stiffness k; by means of which the system can be 
forced when required. Suppose that the angular displacements of the rotors 
are 0,, 62, etc. 

In the absence of damping and with no restraint imposed upon the exten- 
sion shaft, the free motion of the system is given by a set of five equations of 
which the following is typical: 


10 + K(O2 — 93) + k(O2 — 6) = 0 (16.1) 


The system has four principal modes. For any of the natural frequencies 
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b/2m, we may substitute —b6, for 6,, —b?0, for 9,, etc. to obtain a further 
five equations of the form 


—k0, + (2k — Ib?)0, — k0, = 0 (16.2) 
Equations (16.2) may be solved for b in a number of ways, all of which 
are tedious, to yield the following non-trivial values: 
b, = (0:099 k/1)*, bs = (2:18 k/D)* 
by = (0°813 k/D?, bs = (3-46 k/D} 
These give the first, second, third and fourth mode frequencies respectively. 
The relationship between the angular amplitudes of the rotors depends 


on which mode is excited. Equations (16.2) are satisfied by the values given in 
Table 16.1, where the displacement at rotor 1 is taken as unity and positive. 


(16.3) 











TABLE 16.1 
Displacement Mode 
at rotor no. 
first second third fourth 
1 1-000 1-000 1-000 1-000 
2 0-904 0-187 —1:056 —2-439 
3 0-707 —0-778 —0-810 2-562 
4 0-457  —1:102 200 oe 
5 —0-155 0-035 —0-014 0-010 


nnn 


A negative sign indicates where a rotor moves in anti-phase with rotor 1. 
Where a sign change occurs between adjacent rotors, this indicates a node in 
the connecting shaft; there are as many nodes as the order of the mode. 
We may draw the elastic curves as shown in Figure 16.2, using the values 
given in Table 16.1. The ordinate is angular displacement in any convenient 
unit. Straight lines are drawn between the plotted values since each section 
of shaft is assumed to twist uniformly along its length. The nodes occur 
where the lines cross the horizontal axis, the actual positions being found 
by simple proportion. 

For the experiment, the same elastic material with the same cross-section 
should be used for each portion of shaft. The section between rotors 4 
and 5 will be twice as long as each of the other sections in order that its 
stiffness shall be half as much. Hence, in Figure 16.2, the slope of each line 
is a direct measure of the maximum shear stress occurring in the section. 
In practice we must guard against the danger of the effect of resonance, 
at any of the mode frequencies, permitting excessive shear stress and thus 
failure of the shaft. 

Now suppose that a sinusoidal displacement A sin wt is imposed at the 
end of the extension shaft. Theoretically, in the steady state, each rotor 
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will oscillate at the forcing frequency w/27 with some amplitude dependent 
upon A, o, I, k, k, (and, in general, the rotors will have different amplitudes 
from each other). Each rotor will move in phase or in anti-phase with the 
forcing function. When w = by, bg, etc., rotor 1 will have amplitude A and 


S 


Angular displacement of shaft 





FicureE 16.2 Elastic curves for multi-mass system. 


there will be no stress in the extension shaft; the other rotors will move 
according to the relative values in Table 16.1. Effectively the system will be 
vibrating freely. 

Of course, in practice there will be some frictional damping which will 
require a small power input to maintain any of the free modes. This will be 
provided by the source of the forcing function. The effect of the inherent 
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FiGureE 16.3 Forced oscillation at a rotor. 


damping on the natural frequency values is likely to be small. The amplitude 
at each of the rotors may be expected to vary with w in the manner shown in 
Figure 16.3, rising to a peak at b,, bo, etc. 

Were we to impose a sinusoidally-varying torque instead of a displacement, 
the amplitude at resonance would increase theoretically to infinity. There is 
usually insufficient damping in practice to avoid substantial oscillations of the 
shaft at the critical frequencies, followed by shaft failure in extreme cases. 


~ 
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Experimental Details 


For experimental purposes it is desirable to keep the natural frequencies as 
small as possible; this dictates a slender shaft which must be supported in a 
roller bearing close to each rotor. The general arrangement is shown dia- 
grammatically in Figure 16.4. A steel shaft of diameter 3mm and length 
30 cm has a stiffness k of about 2:2 N m/rad. This is suitable for the sections 
between rotors 1 and 2, 2 and 3, 3 and 4, whilst that between 4 and 5 is 
made 60 cm long. The largest rotor is of steel, of diameter 23 cm and thick- 
ness 2cm; this has a mass of about 6:5 kg and a moment of inertia 20 J 
of about 0:042 kg m?. Corresponding smaller rotors each have a moment of 
inertia J of about 0-0021 kg m?. By equation (16.3), the natural frequencies 
should be about 1-6, 4:5, 7-5 and 9-4 Hz. 





Ficure 16.4 Arrangement of experimental system. 


The other advantage of a slender shaft is that large amplitudes may be 
excited without overstressing the shaft. Even so, for the suggested shaft, 
the total twist between any pair of smaller rotors should be limited to 13° 
or 14°. Light pointers may be attached to each rotor to sweep over fixed 
scales marked over a range of +8°. The pointers can be illuminated by a 
stroboscope set to flash at near the frequency of oscillation, so that the pointers 
appear to oscillate at the beat frequency, which can be made small enough to 
enable the rotor amplitudes to be read from the scales. It will also be evident 
whether rotors are oscillating in phase or anti-phase. 

An arrangement is required to turn the end of the extension shaft with 
simple harmonic motion with a maximum amplitude of about +6°. One 
possibility is shown in Figure 16.5, a lever of about 10 cm in length being 
attached to the shaft and moved via a long connecting rod by a crank driven 
by an electric motor. A simple means is required for adjusting the crank 
throw from 2 mm to 1 cm. Though a 0-5 kW variable-speed d.c. motor will 
suffice, a Velodyne (speed-controlled) motor is preferred to minimize cyclic 
speed fluctuations. Accurate speed measurement is essential and may be pro- 
vided by a digital counter driven by a magnetic or a photo-elastic pick-up 
on the motor shaft. 

If desired, the surface stress fluctuation at one or more of the shaft sec- 
tions may be detected by a rosette of strain gauges cemented to the shaft 
with a suitable amplitude-modulated measuring system. 
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It is instructive first to find the theoretical natural frequencies from 
equation (16.3) and to determine the corresponding motor speeds. 

With the crank throw set at minimum, the motor speed may be slowly 
increased until the first mode is excited. With the aid of the stroboscope, 
the amplitudes of the rotors may be checked. Repeat for the other modes in 






Motor shaft 


j Test shaft 
FicurE 16.5 Means to excite the system. 


succession. Use this information to determine the maximum crank throw 
which may be used at each mode to give the largest possible amplitudes with- 
out overstressing the shaft anywhere, and use these values in the subsequent 
tests. 

Concentrate your attention on frequency ranges in the vicinity of the 
four actual critical frequencies. At each setting of the speed control, measure 


Twist /length 





Frequency 


FicureE 16.6 Twist per length in a section of shaft. 


the speed and all the rotor amplitudes. For each shaft section, determine the 
twist per unit length and plot against frequency in the manner of Figure 16.6. 
Compare the peaks in your curves with the theoretical natural frequencies. 
Where strain-gauge equipment is provided, measure the stress and plot 


against frequency. Compare with the relevant twist-frequency plots and 
check by simple torsion theory. 
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Discussion 


Make sure that you can obtain the expressions for the natural frequencies 
given in equation (16.3). State which method you prefer and why. Check 
the values given in Table 16.1. 
Discuss fully the experimental curves. How do the actual critical frequencies 
-differ from the theoretical values? Give reasons for any substantial 
discrepancies. 


7 ics teina 
Speed Control System 


Our last control experiment is concerned with the control of the speed of a 
shaft driven by a prime mover. It is instructive to compare and contrast its 
behaviour with that of the position control system considered earlier. 


Theory 


Figure 14.1 used for the position control system also represents the present 
process. The important difference is that 6, now means the control speed 
of the driven shaft instead of the controlled position. J is the moment of 
inertia of the rotating parts referred to the driven shaft and c is the coefficient 
of the viscous damping. The prime mover in the suggested experiment is a 
reciprocating internal-combustion engine and J and c may reasonably repre- 
sent its dynamics: L is the external load torque to which the prime mover is 
subject. 

Likewise the block diagram of Figure 14.2 can be used to represent the 
closed-loop control system. 6; is now the desired speed. We may write the 
following equations [1]: 





F 75) —L=(D + c)6, (17.1) 


6=0,—9, (17.2) 


Note the difference in the right-hand sides of (17.1) and (14.1) since 8, is 
now a speed. The operator G, is the coefficient of the speed error on the 
left-hand side of (17.1). The proportional control, of coefficient k, is subject 
to the engine induction lag which is approximated to a simple exponential 
lag of time constant 7,, the latter depending on the throttle restriction and 
the manifold capacitance. 

Eliminating 6, from (17.1) and (17.2), we obtain 


{7.[D? + (cr, + DD + (c + AD} 
= {7,JD? + (cr, + DD + c}O; + (l+7,.D)Z (17.3) 


which gives the error as a function of 0, and L. By comparison with equation 
(14.3) for the position control system, (17.3) shows that once again we have 
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a second-order system. In practice the parameters J, c, k and 7, are such 
that the system mode is underdamped oscillatory. The steady-state error, for 
constant L and 6;, is 


O55 = {el(c + kK)}8;, + tlle + kL (17.4) 


.Compared with (14.4), we note that a sustained change in 6, gives rise to 
a steady-state error as well as that given by a load change. 

To avoid steady-state error, it is customary to employ integral control 
action; this means that the controller integrates the error signal with respect 
to time and the new controller output is added to the proportional action. 
Instead of (17.1), we now write 


6 —L=(ID +0), (17.5) 


where 1/D means integration with respect to time and 7; is the integral 
action time constant. It is usually possible to adjust k and 7, separately by 
means of knobs on the controller. 

We now eliminate 0, between (17.5) and (17.2), to obtain 


{7,JD® + (cr, + I)D? + (Cc + AD + k/736 
= {7,JD? + (cr, + DD? + cD}0,+ (1 +7.D)DL (17.6) 


The system is now of third order, there being an aperiodic mode in addition 
to the principal oscillatory mode. The steady-state error, for constant L 
and 6@;, is now 


6,, = 0 (17.7) 


The advantage of integral control action is self-evident. The disadvantage 
is that the system can easily become unstable if k is too great or 7, too small. 
It might be argued that good stability can be restored if a third control action, 
namely derivative or error, is added. Although of great benefit in many 
systems, it is useless in the present case for the following reason. There is an 
inherent cyclic speed fluctuation in a reciprocating engine, the frequency of 
which in a 4-stroke 4-cylinder design can be easily as high as 150 Hz compared 
with the principal system mode frequency of about 0-5 Hz. Differentiation 
of the high-frequency ripple in the error signal amplifies its magnitude by 
2a x frequency. This results in excessive ‘noise’ associated with the controller 
output signal and control is lost. 


Experimental Details 


Test Equipment. It is desirable to avoid non-linearities as far as possible. 
For this reason, the common inertia ‘governor’,* which inherently provides 
an output proportional to the square of speed, should not be employed. 


* In control terminology, the term ‘governor’ means a complete closed-loop speed- 
control system. The ‘governor’ referred to in texts on machines refers only to the speed- 
measuring element of the system. 
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Instead a tacho-generator should be driven by the engine shaft. This has an 
output in the form of a d.c. electric potential which is proportional to speed 
6,; a suitable type gives 0-04 V/(rev/min). This is fed to a pot-coil relay 
to give a current signal of 0-05 mA/V input. 

It is convenient to employ a commercial proportional-plus-integral elec- 
tronic process controller with possibly a small modification to extend the 
7, range down to 1 second. This includes the system difference element, and 
hence the pot-coil relay output is fed in directly, zero to 5000 rev/min in 0, 
changing the controller input over its entire range of zero to 10 mA, moving a 
pointer over the input scale marked 0-100. A knob on the controller varies 
the desired speed 6; fed in, there being a second scale marked 0-100. The 
speed error 0 is shown by the difference between the two scale readings. 

The controller output is also in the range zero to 10 mA and is shown on a 
third scale marked 0-100. The proportional action is usually calibrated as 
proportional band, which means the percentage of the maximum input 
change required to cause full-scale output change. Thus, if a difference in the 
desired and measured value pointers of 50 results in an output scale change 
of 100, the band is said to be 50%. The band may be calibrated from 25 to 
300%; the parameter k of the system is, of course, inversely proportional to 
the band. 

On proportional action alone, it is essential that the controler output 
should be at about mid-range when the desired and measured pointers 
are coincident; this permits changes in output for both positive and negative 
speed errors. The controller should be checked and the zero-error output 
value noted. This is of no consequence when integral action is used. 

A four-cylinder spark-ignition engine, of 1200 cm® displacement, has been 
found suitable for this experiment. The external load would normally be 
provided by a hydraulic dynamometer but, for the purpose of the experiment, 
a means is required for rapidly turning the load-setting handwheel. 

The throttle valve of the single carburettor is required to be turned in 
direct proportion to the controller output. It is convenient to use a small 
pneumatic motor fed via an electro-pneumatic converter, both of which need 
clean air supplied at about 130 kKN/m? gauge pressure. The lag introduced by 
this regulating unit is small enough (compared with the main time constants 
of the system, I/c and 7,) to have negligible effect on the system behaviour. 
Its overall characteristic is linear, as shown by Figure 17.1. 

The equipment is completed by installing a two-pen strip-chart current 
recorder, one pen recording the same as the pointer on the controller input 0, 
scale and the other that of the controller output scale. The record will then 
indicate the variation in the controlled speed together with that of the throttle 
opening, both against time. 

It is also convenient to have a manual control knob on the controller, 
together with a manual/automatic changeover switch. This permits the control 
loop to be broken and the engine throttle to be adjusted directly by hand. 


Engine and Dynamometer Characteristics. The only serious non-linearities 
in the experimental system are due to the engine and dynamometer. Typical 
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Ficure 17.1 Regulating unit characteristic, 
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Ficure 17.2 Engine and dynamometer characteristics. 
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curves are shown in Figure 17.2. The full curves indicate the load imposed by 
the dynamometer for various fixed settings of the handwheel (measured in 
handwheel turns). The dotted curves show the torque output of the engine at 
various fixed throttle positions. At steady speed, the engine output equals the 
dynamometer load. . 

It is evident that the engine output torque is not proportional to throttle 
Opening at any speed. In fact, the system parameter k, though nominally 
fixed by the proportional band of the controller, varies both with speed and 
throttle opening. In order to ensure stable operation under all conditions, 
the band must be set so as to provide satisfactory response at the speed and 
load where the spacing of the throttle curves is greatest. The largest variation 
occurs at low speed, k being much greater at low load than at high load. 

The only convenient way to impose a change in load L is by re-setting the 
dynamometer handwheel. The load set is constant only in the absence of 
a change of speed. Hence, during a transient following a step load change, 
we must realize that the actual load L as well as the speed will fluctuate. 


Procedure. First allow the engine to warm up with moderate load on 
manual control. Set the controller to work on proportional action only with a 
proportional band between 100% and 150 %. Set the dynamometer handwheel 
at about 5 turns. Set the desired value at about 2000 rev/min and switch to 
automatic control. If the system tends to hunt, switch back at once to manual 
control, raise the proportional band to stabilize and repeat. Now impose 
small changes in 0; and check whether the response is oscillatory but reason- 
ably damped. Otherwise adjust the band further. The system is now ready 
for the following suggested tests. 


(a) With the dynamometer handwheel set at 5 turns and the proportional 
band at the value chosen above, set 0; to 1500 rev/min (30 on the scale) 
and allow the system to settle. Start the recorder. Make a step change in 
0, to 3500 rev/min (70 on the scale). Allow 10 seconds for the system to 
settle completely. Now re-set 6, to 1500 rev/min and repeat. You should 
now have a record on the lines of Figure 17.3. 

Repeat for different loads and steps in 0;. For high load and large 
steps, note that the system may saturate during the early part of each 
transient, i.e. the throttle may be held at a maximum or minimum 
opening. 

(b) Re-set the load at 5 turns. Add integral control action of time constant 
of about 1 second and re-adjust the proportional band if necessary. 
To avoid saturation, use smaller steps in 6; than in (a). For example, 
start at 1800 rev/min and make a step change to 3200 rev/min. Figure 
17.4 shows the kind of response obtainable. Note the removal of 
steady-state error. Repeat for different loads. 

(c) With the same controller settings as the above, set 0; at about 2800 
rev/min and the dynamometer handwheel to 1 turn. After the system 
has settled, apply a step change in load by turning the handwheel 
rapidly to 6 turns. Allow 10 seconds for the system to settle and then 
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Ficure 17.3 Response to desired-value change with proportional control. 
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FicurE 17.4 Response to desired-value change with proportional-plus-integral 
control. 
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Ficure 17.5 Response to load change with proportional-plus-integral control. 
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apply a reverse step in load back to 1 turn. A typical response is shown 
in Figure 17.5. Repeat for different speeds and steps in load. 


Discussion 7 


Check that you can obtain equations (17.3), (17.4), (17.6) and (17.7). Draw 
a detailed block diagram of the experimental system, identifying all the 
components employed. 

Explain the large discrepancies between desired and actual speed in the 
steady state in test (a). In test (b), following the step up, why is there a 
well in the controller output response followed by a gradual rise to the final 
value? 

Why is the response for the step down more oscillatory than that for the 
step up in all three tests? 


Reference 
1. WEBB, C. R., Automatic Control, McGraw-Hill, 1964. 
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Analogue Computation 


An analogue computer is an instrument used for making a mathematical 
model of a system. Unlike a digital computer, which is fast and accurate and 
works in small discrete steps, an analogue computer represents variables 
continuously. It solves the differential equations representing the system 
directly, and hence is easier to program than a digital machine. It shows to 
advantage when simulating complex systems which are laborious to analyse 
by ordinary methods. Its accuracy depends on that of its components and on 
the skill of the programmer. 

We shall confine our attention to the kind of instrument which is based on 
the d.c. operational amplifier. 


Computer Elements 


When the output of a high-gain d.c. amplifier is fed back to the input, it acts 
as an amplifier of gain —1 if the input and feedback resistances are equal. 
If the former is smaller than the latter, the magnitude of the gain is greater 
than unity in inverse proportion. Variants of the above are used for addition, 
subtraction, integration with respect to time, etc. Table 18.1 shows three 
important forms with their usual symbols. The first can add a number of 
voltages, multiplying inputs as required by round factors such as 10 or 100; 
the sign reversal is inherent. The second is a simplified special-purpose form 
of the first. The third element integrates with respect to time and permits 
initial conditions (I.C.) to be inserted. 

Another important element is the attenuator, which is also shown in 
Table 18.1. This is simply a precision potentiometer by which the output 
voltage can be made some fraction / of the input. 


A Specific Problem 


The use of a computer is best described with the aid of a particular example. 
We choose a linear mechanical vibrating system which is just sufficiently 
difficult to analyse to justify the use of a computer. The same system is used 
in the actual experiment on the computer. Typical commercial instruments 
can, of course, deal with much more complex systems which can be non- 
linear or be represented by partial differential equations. 
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TABLE 18.1 





Symbol Function 





V4 Summing amplifier 


fe) Vo = —(Y, + 10V2) 


Vo 


Inverter 
V, Vv 
; a Vo=—-V, 
Y 
Integrator 


2 - 


Attenuator 


Vi Vo=4Vv, (A=1) 


~ 


be 


Vo 


Figure 18.1 shows a lumped system consisting of two bodies, two springs 
and a viscous dashpot. The following data apply: 


m, = 25 kg 
Mz = 2kg 

k, = 10 N/cem 
kg = 50 N/em 


c = 1 N/(cm/s) 






kp 





iy, it gig obo 


FicureE 18.1 Lumped mechanical system. 
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The large differences in the chosen parameters are deliberate in order to 
illustrate the versatility of the computer. We are interested in the displace- 
ments x, and x, resulting from an input displacement x3. The equations of 
motion may be written as under 


mX, + c(% — Xe) + kiG — x2) = 0 (18.1) 
MoX_ + C(X_ — X1) + ka(Xa — 1) + ko(x2 — Xs) = 0 (18.2) 


The first step is to re-write the equations with the higher derivatives on the 
left and the remainder on the right, i.e. 


XH, = —(c/mi)x% + (clm)X2 — (kal). + (Ki /11) x2 (18.3) 


Hq = —(c/me)X2 + (c/ma)X, — (ky/im2)x2 
+ (ky/ms)x1 — (K2/me)x2 + (K2/mz)x3 (18.4) 


Amplitude Scaling 


We must not exceed the maximum voltages permitted in the computer; 
otherwise certain circuits will overload (indicated by warning lights) and the 
problem will not be solved. Assuming that we may work in the range zero to 
10 V, positive and negative, we now have to estimate the maximum likely 
values of the variables and their derivatives. 

It is proposed to investigate the transient response in x, and x, to a step 
in x3 of 1 cm. It is then likely, for example, that x, will not exceed 2 cm 
and hence we may represent 2 cm by 10 V. We therefore apply a scale factor 
of 5 V/cm and represent x, by the computer variable [5x1], which of course is 
measured in volts. Other computer variables are found likewise, as in Table 
18.2. 





TABLE 18.2 
Problem Maximum Scale Computer 
variable value factor variable 
eee ene a Sine Se 

Xe 2cm 10/2 [5 x1] 

Xe 2em 10/2 [S x2] 

Xe 1cm 10/1 [10 x3] 

X41 5 cm/s 10/5 [2 Xi] 

Xe 25 cm/s 10/25 [0-4 x2] 


oe ie 


Since %, and X, are not to be found, they may be formed as integrator 
inputs. They may be scaled as we please but it is convenient to use the same 
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scales as for x; and X2 respentivelyy giving the computer variables [2*,] and 
[0-4%.]. 

We now re-write the equations of motion in terms of the computer variables, 
1.€. 


S 


pe] = — = pal + 5~ fot — 5 bm] + 3 “1 15x5] (18.5) 


rae = pe] — 3 [Ba 


2 ky 


7 50 


[5x3] — a [Sx] +55 “2 [10x (18.6) 


Time Scaling 


Using the data, we may find the coefficients of the computer variables given 
in the above equations. These vary in size from (c/m,) which is 4 units to 
#5 (K2/mz) which is 200 units. We shall assume that the operational amplifiers 
in the available computer provide gains of either 1 or 10. Clearly it will reduce 
the number of computer elements needed if we can reduce the magnitude of 
the larger coefficients. This will automatically reduce the speed at which the 
computer solves the problems. The process is known as time scaling; we 
apply a factor a to all the coefficients in (18.5) and (18.6). It should be noted 
that, for a process which is slow in real time, a fractional value can be chosen 
for a in order to make ‘computer time’ faster than real time. For the present 


TABLE 18.3 


Number in 





Attenuator equations Expression A 
number (18.7) and represented value 
(18.8) 

1 (i) c/(am,) 0-40 

Z (ii) c/(2am,) 0-20 

3,4 (iii), (iv) ky |(25amy) 0:16 

5 (v) c/(10amz) 0-50 

6 (vi) c/(Samz) 1-00 

7,8 (vii), (viii) k,/G25amz) 0-40 

9 (ix) ke/(1250amz) 0-20 


10 (x) k2/(250amz) 1-00 
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problem, a suitable value for a is 10. The equations are once again re-written 
as 


1 
= [2x1] = —1 (i) [241] + 10 (ii) [0-4%2] — 10 (iii) [5x] + 10 Giv) [5x0] 
(18.7) 


1 
-, [0-42] = —10 (v) [0-42] + 1 (vi) [241] — 10 (vii) [5x2] 


+ 10 (viii) [5x] — 10 x 10 (ix) [5x2] + 10 (x) [10x3] (18.8) 


In the above, the unbracketed coefficients represent amplifier gains, 
either | or 10, The terms in curved brackets are attenuator settings 4; these 
must not exceed unity and dictate the choices of amplifier gain. The values of 
A are listed in Table 18.3. 

We may now draw the computer block diagram, Figure 18.2. Integrators 
1 and 3 have inputs corresponding to the right-hand sides of equations 
(18.7) and (18.8) respectively. Hence the input to integrator 1, for example, 
represents 2%,/a, i.e. 


: 2d 
aoe 
inpu adt (X1) 
ee 
=25 @) 
where 7 is computer time, equal to at. Therefore 


output = —[2x,] 


which is the negative of the required computer variable. By a similar argu- 
ment, the output of integrator 3 is [0-4x2], all the terms on the right-hand 
side of (18.8) being fed in negatively for convenience. 

The output from integrator 2 is required to be [5x,], and hence the input 
must be —5x,/a, i.e. —4x,. We must therefore insert attenuator 11 of value 
0-25. Similarly attenuator 12 in the input to integrator 4 must have a value of 
0-125, together with an integrator gain of 10. 

Inverters 5, 6 and 7 are needed to reverse signs in certain paths. At this 
point it should be stressed that the arrangement of Figure 18.2 is not the 
only one possible. However, it is important to aim at minimizing the number 
of computer elements required. 


Static Check 


The object of this part of the procedure is to check the programming. We 
make use of the facility of applying initial conditions (I.C.) to the integrators. 
These need not be realistic but are chosen to ensure that voltages which are 
neither zero nor greater than 10 V in magnitude are produced in every path 
in the diagram. We assume the values given in Table 18.4. 
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96 





Ficure 18.2 Computer block diagram. 


ANALOGUE COMPUTATION 97 
TABLE 18.4 


Variable Value of Corresponding value 


LG of computer variable 
X1 lcm 5 V 
Xe —lcm —5V 
xy 5 cm/s 10 V 


Xe —10 cm/s —4V 


For the purpose of the check, only the I.C. signals operate the integrators, 
the other inputs being switched out. We require —4 V at the output from 
integrator 3, for example, and must therefore apply 4 V; this is done by 
applying 10 V (marked ‘S.C.’) through attenuator 13 having value 0-40. 
Similar S.C. signals are applied to integrators 1, 2, 4. Another 10 V signal 
is used to represent the input [10x3] (marked “Input’). 

We now trace through the computer block diagram, checking that a 
voltage other than zero but not exceeding 10 V appears at every point. If 
these requirements are not met, the assumed initial conditions must be 
changed accordingly. The aggregate inputs to integrators 1 and 3 are then 
found to be —28 V and —270 V respectively. These are checked by inserting 
the assumed initial conditions in the original equations (18.3) and (18.4). 
For example, 


¥, = —4(5) + 4(—10) — 40(1) + 40(—1) 
= —140 cm/s? 


1 
Hence = [2x,] = —28 V, as required. 


The above procedure checks the whole of the programming and suggests 
1.C. values suitable for an actual static check on the computer itself. 


Setting Up the Computer 


This varies in detail from one instrument to another and reference should 
be made to the manufacturer’s literature. However, the general procedure 
is as follows. First make the necessary connections on the patch board to 
represent the computer diagram, Figure 18.2. Next adjust all the potentio- 
meters to their required values. 

Now carry out a static check on the same lines as described above; in 
most instruments special switching arrangements are available for this. The 
instrument should not overload anywhere and all the required voltages 
should be obtained, provided the setting up has been carried out correctly. 

Make connections to the recording equipment to record the variables 
[5x3], [5x2] and [10x]. 
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Computer Experiment 


Switch out the initial conditions used for the static check and switch in the 
input variables to the integrators. No initial conditions are required for the 
test considered here, though of course transients can be started from finite 
initial conditions whenever this is required. 

Use the appropriate controls to impose a 10 V step in [10x3] which is 


Magnitude ratio 





t (s) 


FiGureE 18.3 Step response of mechanical system. 


equivalent to.a 1cm step in x3. The computer should then solve for the 
transient response in x, and x2, producing records of [5x,] and [5x,] ten 
times slower than real time. These records should be translated into actual 
variables and real time so that the ratios |x,/x3| and |x2/xs| appear as in 
Figure 18.3. 

Further tests can now be made to find the effect of changing the damping 
coefficient c. For ¢ equal to 0-5 N/(cm/s), for example, it is only necessary 
to halve the values of attenuators 1, 2, 5 and 6. The instrument may overload, 
however, should the changes cause maximum transient values exceeding 
those assumed in setting up the program. : 


~ 
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Discussion 


Make sure that you understand the functions of amplitude and time scaling. 
What is the advantage of slowing down the solution to the given problem? 
Try to find out exactly what switching operations are performed by the 
controls provided on the available instrument. How accurately can you 
- set the attenuator values? What lack of accuracy was revealed by the static 
check on the instrument? Are you satisfied that your final curves represent 
an acceptable solution of the mathematical problem? 
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Elastic Constants 





Before any useful calculations or experiments on stressed bodies may be 
performed, the stress-strain relation of the materials must be determined. 
Experimentally, we can measure only strain or deflection, from which the 
stresses may then be inferred. 

For most common structural materials, the stress-strain relation is linear, 
or almost so, in the elastic region. This region is characterized by no per- 
manent deformation after load is removed, i.e. complete return to original 
size. Two types of stress may be applied, normal or shear, and two constants at 
least must be determined. This experiment will determine one constant, and 
another will be determined in the succeeding experiment. 


Theory 
There are four elastic constants in common use: 
(a) Young’s modulus, E, defined by 
normal stress 
normal strain 
(b) Rigidity, or shear, modulus, G, defined by 
shear stress 
shear strain 
(c) Poisson’s ratio, v, defined by 
lateral strain 
normal strain 
(d) Bulk modulus, K, defined by 


hydrostatic stress 
volume strain 


The first three are all used in structural calculations. If the material under 
test is isotropic, i.e. its properties are the same in all directions, and homo- 
geneous, i.e. its properties are the same throughout its volume, the four 
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constants will specify the behaviour completely. Two independent equations 
may be found relating them: 


E = 3K(1 — 2p), E=2G(1 + ») 


so that an isotropic, homogeneous structural material is completely specified 
by two elastic constants. 


Experimental Details for Young’s Modulus 


If a uniform bar, cross sectional area A and length L is extended by 
amount S due to axial load P, then 





E = PL|AS 
Hence S = PL/AE 
Loaded wire 
Unloaded ‘i 
Wits Area of 
uniform 
stress 
1 Gauge points 
Vernier = for extensometer 
Weights 
(a) (b) 


FicureE 19.1 Measurement of Young’s modulus (a) vernier method for measuring 
Young’s modulus of wire (b) tensile apparatus for measuring Young’s modulus 
of bar. 


To make S as large as possible for accurate measurement, L must be 
large and A small. This is the basis of the vernier and Searle’s method for 
determining Young’s modulus, where the specimen is in the form of a 
long, thin wire, up to 5 m long. Variation of temperature along such a length 
would cause erroneous results due to a considerable thermal expansion, and 
so an identical wire is hung close to the specimen wire. To the ends of 
these wires is attached a vernier measuring device, as in Figure 19.1, and on 
application of the load the relative extension may be accurately measured. 

If the material cannot be obtained in wire form, then measurement must be 
made on a machined bar. Again, the bar should be kept as long and slender 
as possible. The load will now have to be applied by a lever or hydraulic 
machine to obtain the same stress levels. Measurement of the extension is 


~ 
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not so simple, as the stress is not uniform at the ends of the bar where it is 
attached to the loading machine. In this case, the extension must be measured 
over a central portion of the bar by an ‘extensometer’. This is a compact 
mechanical gauge which, when clamped to the specimen, magnifies the exten- 
sion of a fixed length known as the ‘gauge length’. Common types of extenso- 
meter are the Huggenberger and Demec gauges, which are available with 
different gauge lengths. The extension should be measured on opposite 
sides of the bar to provide an average value. 

Plot load against extension and derive a value of Young’s modulus for 
brass, copper, aluminium and mild steel. Repeat these measurements for the 
materials in bar form. 


Discussion 


Compare the results for materials in wire and bar form and comment on 
them. Discuss the sources of error in both experiments and suggest how 
they may be reduced: estimate the accuracy of the results. 
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Tapia of a Circular Shaft 





Torsion is acommon type of load applied to structural members. When applied 
to non-circular cross-sections, the analysis is complex, so this experiment is 
restricted to circular shafts. The object will be to verify the torsion formulae, 
and to determine the rigidity modulus for several materials. 


Fixed end 
- 
/ ~ 





Rotated end 


Ficure 20.1 Torque acting on a circular shaft. 


Theory 


For a circular shaft under pure torque T, illustrated in Figure 20.1, we make 
the following assumptions: 


(i) twisting is uniform along lengths of shaft, 
(ii) radii remain straight, and 
(iii) cross-sections remain plane. 


The first two assumptions follow if the shaft is elastic, while the third is a 
consequence of the circular symmetry. 
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Considering an elemental tube, radius r, thickness ér, inside the shaft, 
we see that 


shear strain y = r6/L 
But 
rigidity modulus G = z/y 
Hence 
+ = Gr6/L 
Now assuming constant shear stress 7 in the elemental tube, the torque on 


the tube is given by (27rér)rr = (27 ror) = r 


G0 
= 2n— rr 


L 


Hence the torque 7 on the whole shaft is given by 
a GO 
i = 
T= { 2a rae dr 
G6 [r*]* 
aac 
ee E 
Tye) Sey 


The expression in brackets is the ‘polar moment of inertia of the cross- 
section’, usually represented by J. 


Experimental Details 


The circular shaft should be machined as in Figure 20.2. One end is rigidly 
clamped, and a radius arm is fixed to the other, so that torque may be 
applied by weights. To prevent any flexure, the shaft should be supported 
by a tail-stock at its free end. 

Because of the non-uniform cross-section, the angular deflection is best 
measured by a torsion meter clamped to the middle of the shaft. These 
meters consist of two clamping rings a fixed distance apart (the gauge 
length), which measure the deflection by the relative movement of two 
levers. The length L in the equation is now the gauge length of the meter. 

To allow for any initial slip in clamped parts, it is advisible to apply a 
small load before measurements are made. As an extra check for this type of 
error, readings should be taken for loading and unloading. 

Plot load against deflection and derive a value for the rigidity modulus. 
Do this for specimens of brass, copper, aluminium and mild steel. 
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Clamped end 





( = 25cm 


mtn 


Hole for locating pin (tail stock) 
25cm 


FiGurE 20.2 Specimen shape for torsion experiment. 


Discussion 


Discuss the sources of error in the experiments, and estimate the accuracy 
of the results. Describe what would happen if the elastic limit were exceeded, 
and consider how well the three assumptions made in the theory have been 


justified. 


Combine the results with the appropriate Young’s modulus (obtained from 
Experiment 19) and plot E against G. What conclusions can you draw from 


this graph? 


21 
Simple Bending 





Bending of a beam is nearly always accompanied by shear forces, so that in 
addition to the bending stresses we have shear stresses. There is one type of 
loading (commonly used in testing) known as four-point loading, which 
produces a region of pure bending in the beam. This enables us to verify 
the equation for pure bending and determine the second moments of area 
experimentally. 


Y Area 5A 






Tension 





Compression 


CROSS-SECTION STRESS DISTRIBUTION 


Ficure 21.1 Element of a beam under the action of a bending moment M. 


Theory 
Consider a small element of beam under a pure bending moment, as in Figure 
21.1. We will assume: 

(i) Young’s modulus is the same in tension and compression, ' 

(ii) the element was initially straight, and all filaments bend into circular 
arcs with a common centre of curvature, the radius of curvature R 
being large, and 

(iii) transverse cross-sections remain plane. 
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The first assumption is true enough for most elastic materials, and the 
second is, geometrically, a good approximation for large radii of curvature. 
The third assumption is not true if shear forces are present, as these distort 
the transverse planes; the distortion is small, however, and this derivation 
may be used when shear occurs. The cross-sectional shape is also distorted 
by the Poisson strains, caused by the varying bending stresses, and known 
as anticlastic curvature. This does not affect our analysis. 

Assuming a neutral axis in the element, where the strain and stress are zero, 
then strain in a filament distance y from this axis is given by 

(R+y)0—RO o y 
RO "ie ak 
From the equilibrium of forces, we have three conditions: 
(a) A small filament, cross-section 6A, distance y from the neutral axis 
will be in equilibrium so the nett normal force on the total cross- 
section is zero, i.e. 


fodA =0 
E 
or R [> dA =0 
Thus the neutral axis must pass through the centroid of the cross- 


section. 


(b) The bending moment is balanced by the moment of normal forces 
about the neutral axis, i.e. 


M = foydA 
E 
_ — 2 
=R | y?dAa 
= El/R 
where J is the second moment of area. For a rectangular beam 


1 
I= = x breadth x (height)°. 


12 
(c) For equilibrium about axis YY (perpendicular to the neutral axis) 
foxdéA =0 
E 
or R | yx ddA =0 


This integral is the product of inertia, and the axes for which it is zero 
are the Principal Axes of Bending. These axes are parallel to the sides of 
a rectangular cross-section. This point will be exemplified in a later 
experiment. 


We thus have the equations of bending: 
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Experimental Details 


The particular loading arrangement is shown in Figure 21.2, along with the 
shear force and bending moment diagrams. The bending moment is constant 
between the supports, so the radius of curvature is also constant over this 







Supports 


Weignt W 


Shear force 


Bending 
moment 
wa 
2 





Ficure 21.2 Loading scheme, with associated shear-force and bending-moment 
diagrams, for a simple beam. 


length. Thus the beam traces out the arc of a circle, illustrated in Figure 
21.3, and this is termed circular bending. 

If 6 is the upward deflection of the centre point B, then from the geometry 
of intersecting chords 


ORS = L? 
1a 
Me REI 
# yee Et 
ence = EI 


If the depth of the beam is 2h, and the maximum strain at section Bis €max, 
then 


Omax = 2z@max = Mh|I 
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The quantity I/h is called the section modulus and is usually represented as Z. 
Therefore 


Cmax = M/ZE 


Measure the deflection and strain at point B, for increasing load, and repeat 
this for beams with different cross-sectional areas. 





/ 


Radius of curvature R 


FicureE 21.3 Deflection of the centre part of a beam under pure bending. 


It is advisable to use only rectangular cross sections, having depths less than 
breadths. 

Plot load against deflection, and load against strain, and derive values of 
the second moments of area. Compare these experimental values with those 
from direct measurements of the cross-section, giving estimates of error in 
each case. 

Choosing a suitable value of load, plot graphs of deflection, 6, against 
J, and strain, @max, against Z. 


Discussion 


Derive the expression for the second moment of area of a rectangular cross- 
section, and comment on the values determined experimentally. Consider how 
well the results verify the equation of bending, and suggest which assumption 
is most likely to be in error. 
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Bending of an L-shaped 
Cantilever 





In the experiment on simple bending, the beam was loaded parallel to a 
principal axis, so the product of inertia was zero. In this experiment, the full 
bending theory is demonstrated by measuring the deflections of a cantilever 
with an L-shaped cross-section under the same load, but at different orienta- 
tions to the principal axes. The positions of these axes are determined 
experimentally and by calculation. 


Theory 


Consider the irregular cross-section in Figure 22.1, with principal axes Gu 
and Gv, where G is the centroid of cross-section. 


Deflection 5S. 





FiGurE 22.1 Principal axes and deflections of an asymmetrical beam (a) cross- 
section of an asymmetrical beam, (b) resolution of the deflections of an asym- 
metrical beam. 


Let 3A be a small element of area, the co-ordinates being (w, v) with respect 
to Gu and Gv and (x, y) with respect to Gx and Gy. Then 


u=xcosa+ysina 


p= ycosa — xsin & 
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The product of inertia is defined as 
ly = fuvdA 
= J(x cos « + y sin «)(y cos a — x sin a) dA 
= sina cos a[fy? dd — fx? dA] + (cos? « — sin? a)fxy dA 
= sin 2a(/,, — [,,) + cos 2a/,, 


For the principal axes 


I. =0 
pd 
Hence tan 2% = isacie 
Similarly it can be shown that 
Tin fv? dA 
= I,, cos? a + Jy, sin? « — I,, sin 2x 
Ly = Ju A 


= I[,, cos? « + J,, sin? « + J,, sin 2« 
uy zx Yy 


Consider a cantilever with end load W. Let W be inclined at an angle ¢ to 
the principal plane Gu. The deflection along Gu is given by 


W cos ¢. [8 
= §ET = O1, say 
and the deflection along Gv is given by 


W sin ¢ . [° 
oe = 06g, say 


where / is the length of cantilever and 
E is the Young’s Modulus of the material. 


The deflection in the direction of load is given by 


WIS 
6, cos¢ + 62 sing = 3EIT. (Luu Cos? 6 + Ipy sin? d) 


The deflection in the direction perpendicular to the load is given by 
‘ WI? sin 24 
6, sind — 6. cos = “3EII,, weit 9; 
Experimental Details 


The L-shaped beam should be approximately 2-5 cm on its side and 1 m long. 
Each end should be attached rigidly to a wheel, one wheel being fixed to a 


as 
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frame capable of rotation against a scale calibrated in degrees. The other 
wheel carries a slide from which the weight is suspended. The whole apparatus 
is illustrated in Figure 22.2. Two dial gauges are fixed against the load 
carrying wheel, in directions corresponding to zero and 90° of the scale. 
The cantilever is rotated by fixed angular increments, say 22:5°, and the 


Pivot 





Encastre end 


FicureE 22.2 Cantilever scheme for an asymmetrical-beam experiment. 


dial gauges read for the addition of a fixed load. It is best to apply the load 
each time, to allow for rigid body movements in the apparatus. © 

Plot graphs of the 0° and 90° deflections against angle of rotation. Deter- 
mine the positions of minimum deflection, and hence determine the values of 
the two second moments of area. Calculate the value of these moments 
from the cross-section dimensions, and compare the two results. 


Discussion 


Suggest what sources of error affect the experimental result, and consider 
what approximation the theory makes with regard to the shear force. 
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Compound Stress and Strain 


The states of stress and strain in a biaxial stress field may be expressed by 
three parameters, usually taken as two perpendicular normal components 
and a shear component of either stress or strain. These may then be used to 
calculate the principal stresses (and strains) and their directions. Provided 
the material is isotropic, the principal stresses and strains have the same 
directions, so that once the principal strains are determined, the stresses 
may be derived very simply using known values of the elastic constants. 

Strain measuring equipment is sensitive only to extensional strain and, 
as shear strain cannot be measured directly, the biaxial state of strain in 
a particular problem must be determined from at least three measurements of 
extensional strain. This problem applies particularly to the surface strains of 
any structure, as the third principal strain is always normal to a free surface. 


Theory 


Analysis of Strain. Suppose e,, ey and y,, are the linear and shear strains in 
the plane Oxy. It is required to find an expression for é,, the linear strain in a 
direction inclined at 0 to Ox, in terms of e,, e,, yz,, and @. 

In Figure 23.1, OP, of length r, is the diagonal of a rectangle, which under 
the given strains distorts into the dotted parallelogram with P moving to P’. 
Remembering that actual strains are very small, 


PP'=PQcos@ + QRsin 6 + RP’ cos 6 
= (rcos 0. e,) cos 6 + (rsin 6. e,) sin 8 + (r sin 8. y,,) cos 0. 
= re, cos? 6 + re, sin? 6 + r,, sin 0 cos 6 


But 
e= PP'/r 
= te,(1 + cos 26) + $e,(1 — cos 20) + 4y,, sin 20 
= 3(e, + ey) + He, — ey) cos 20 + 4y,, sin 20 (23.1) 


This can be compared with the similar equation for stress. The principal 
strains e, and é2, being the maximum and minimum values of strain, occur 
at values of 6 obtained by equating de,/d6 to zero, i.e. 


tan 20 = y,,/(e; — ey) (23.2) 
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Then e, and e, are given by 


(C2 + ey) £ aVi(ez — ey)? + Yiu} (23.3) 


In order to evaluate e,, e,, and y,, (and hence the principal strains) it is 
necessary in general to know the linear strains in any three directions at a 
particular point. If the principal directions are known only two strains are 
required, since 71. = 0 and e, = e, e, = eg. 





FicurE 23.1 Deformation of a rectangular component under the action of three 
strain components, é;, éy and 7zy. 


Stress-Strain Relations. A biaxial state of stress, characterized by principal 
stresses o, and o, will produce the following strains: 


01 02 
o,/E in direction 1 o./E in direction 2 
—vo;,/E in direction 2 —vo,/E in direction | 


where E is Young’s modulus and » is Poisson’s ratio. 
Assuming the strains are small, by the principle of superposition 


strain in direction 1, e, = (0, — v02)/E 
strain in direction 2, e, = (02 — v0,)/E 


Experimental Details 


The thin metal plate illustrated in Figure 23.2 should be loaded by dead 
weights. The particular loading arrangement, known as a ‘whipple tree’, en- 
sures a uniform stress state over the centre of the plate. This area of uniform 
stress should be at least 25 cm square. The plate should have extensometer 
clips fixed parallel and perpendicular to the direction of load, and also at 10° 
intervals between these lines. A Demec gauge with 15 cm gauge length is the 
most convenient instrument for this experiment. To counteract any initial 
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curvature of the plate, extensometer measurements should be made on both 
sides, and the mean reading taken. 

Load the plate in equal increments to just under the elastic limit, and take 
extensometer readings parallel and perpendicular to the load direction. Then 
load the plate somewhere in this range, and take extensometer measurements 
at 10° intervals. 

Plot load-extension curves for the readings parallel and perpendicular to 






Slots and spacers 
all have same width 





25cm 
Extensometer 
pegs 
Ie low. &)4 Rah Oe be a 
| | 
| | 
| ee 25cm SSP east | 
| Clamp | 
et eee eee = =i 


FiGurE 23.2 Plate with slots, for uniform region of strain. 


the load, and determine the elastic constants. Then plot the extensometer 
readings against cos 20, where 0 is the angular position of the extensometer 
location measured from the load direction. As the load direction is a principal 
stress (and hence strain) direction, the shear strain y,, is zero, and the 
equation is 


€o = 3(e1 + €2) + 3(€1 — eg) Cos 20 


From the slope and intercept of this graph, determine e,, e, and hence 
01, O2. Compare experimental and applied stresses. 


- 
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Discussion 


Suggest what sources of error may effect your results and discuss how well 
the experiment justifies the theory. Estimate the size of strains that would 
produce significant errors in the approximations made in the analysis. 

Suggest how the principal strains could be determined if their directions 
were also unknown. 


5—II 


24 
Elastic Stability of Struts 


A thin strut that is compressed along its axis buckles at a certain critical 
load. It is of interest to compare the critical load as obtained experimentally 
with the value predicted by Euler’s theory of thin struts, according to which 
EI 
(— OF 
where P, is the critical load, E is Young’s modulus, J is the geometrical 
moment of inertia of the cross-section of the strut, and L is its length. 
Direct measurement of this Euler load is not easy. In principle, it can be 
measured by adding weights to the top of a vertical strut, and finding the 
load under which any slight lateral deflection given to the strut is maintained, 
as in theory this load is only infinitesimally less than the critical load. How- 
ever, in practice the strut will be neither precisely straight, nor precisely 


vertical, and it is best to determine the Euler load in a somewhat less direct 
manner. 


Theory 


Consider a strut with an initial radius of curvature R, under an end load P 
as in Figure 24.1. From beam theory, 


1 : 
R, = ay, jax approximately 


and the bending moment M at distance x from the top is given by 
M = —Py 
Using the formula for an initially curved beam, 
EK1/R — 1/R,) = M 
then EI, d?y/dx? = M + EI. d?y,/dx? 
or d?y/dx? + ay = d?y,/dx? 
where a’ ae P/EL 
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Initial 
curvature 


Ficure 24.1 Deformation of an initially curved, pin-ended strut, under end- 
loads P. 


The initial shape of the strut y, may be assumed circular, parabolic, or 
sinusoidal without making much difference to the final result, but the most 
convenient form is 


VaR. Sin (orx/L.) 


which satisfies the end conditions and corresponds to a maximum deflection 
c. Any other shape could be analysed into a Fourier series of sine terms. 
Then 


d2y/dx? + ay = —(c. 7?/L”) (sin 7x/L) 
The complete solution is 
he ee ; cm? ]L? . WX 
vy = Asinax + Bcos ax — (PIL) + a ae sin | 


When x = 0, y = 0. Therefore 





(SG) 
When x = L/2, dy/dx = 0. Therefore 
A=0 
Hence 
cm?/L? 1X 
a a? /L? — a? sin L 
2 ere 
apes p sin 


where the Euler critical load P, = 7? EI/L’. 
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The deflection of the centre of the strut is 


cre 
JYmax = Pi—P = A, Say 


Rearranging, 


A S 
A= te P iG: 
so that a plot of A against A/P gives a straight line of slope P,. Such a 
graph is known as a Southwell plot. 

Struts of the same material may be characterized by the term J/L?, which 
may be simplified by expressing J as 


I= Ak? 


where A is the area of cross-section and k is the radius of gyration. 
The term ‘slenderness ratio’ is now defined as 


slenderness ratio = L/k 


and the Euler load can be related to this expression. 


Experimental Details 


The apparatus, due to J. M. Townson, is illustrated in Figure 24.2. A 
simple screw jack is mounted at one end of a base. At the other end is a 
load indicator which consists of a spring with a pointer that moves over a 
scale calibrated in newtons. On the base between these two items is etched a 
lateral deflection scale. The scale is a set of parallel lines at intervals of 1 mm. 
It ranges up to 4cm on each side of the centre line, and is long enough to 
cover the mid-span positions of any strut. Different lengths of struts are 
accommodated by mounting the screw jack at different places on the base, 
sets of holes being bored in the base for this purpose. 

The strut whose Euler load is to be determined must have, or be given, a 
slight curvature in the absence of any load. It has knife-edge ends, and is 
inserted between V-slots in pistons in the screw jack and in the load indicator. 

The best struts to use are of rectangular section, with the depth much 
smaller than the width. Provided that the slenderness and initial curvature of 
the strut are such that a deflection of 2 or 3 cm is obtained at half the Euler 
load, satisfactory plots can be obtained. 

Take rectangular struts of a material with a low Young’s modulus (plastic 
or wood), with slenderness ratios ranging from 100 upwards. Load each 
strut by increments up to about half the Euler load, and measure the maxi- 
mum deflection at each increment. 

Draw Southwell plots for. each strut. Compare the values for the Euler 
load of the struts of equal slenderness ratio, but different Young’s modulus. 
If necessary, determine the Young’s modulus for each strut by a simple four- 
point bending test, and compare calculated Euler loads with those measured. 
Divide the measured loads by the Young’s modulus, and again compare 
these results. 
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Ficure 24.2 Townson’s apparatus, showing load scale, piston, strut and adjustable 
jack. 


Taking the struts of the same material, plot the measured and calculated 
Euler loads against the slenderness ratios. 


Discussion 


Discuss differences between calculated and measured Euler loads for the 
different materials and consider how well the ratios P,/E agree for the 
different materials. 

Explain the significance of the graph P, versus slenderness ratio and discuss 
how well this graph verifies Euler’s theory, accounting for any deviations from 


the Euler theory. 


Resistance Strain Gauges 


The resistance strain gauge is one of the most important developments in 
instrumentation for experimental stress analysis. Invented in the 1930's, it 
quickly replaced mechanical gauges as the standard means of measuring 
strains on actual structures and test models. Though other electrical gauges 
are available, such as the inductive extensometers, the resistance gauge has 
gained popularity because of its cheapness, versatility, and high sensitivity. 
It only provides a measure of the linear (normal) strain at one point, but a 
biaxial strain field may be determined by a combination of three gauges 
about the same point, and whole structures may be investigated using a 
large number of gauges to give a fairly complete stress analysis. Measurement 
can be carried out at one consol. 


Theory 


The electrical resistance strain gauge operates on a principle discovered 
in 1856 by Lord Kelvin, that when an electrical conductor is strained its 
resistance varies in proportion to the strain, extensional strains producing for 
most metals an increase in resistance, compressional strains producing a 
decrease. 

If a conducting wire is bonded to a structure so that strains experienced 
by the structure are transmitted to the wire, a measure of the wire’s change in 
resistance will give an indication of the surface strain in the direction of the 
wire. | 

As the strains experienced by most structural materials in their normal 
working range are very small, the changes of resistance produced in the wire 
will also be small, and will require extremely sensitive instrumentation in 
order to measure them. 

The relationship between deformation and change of resistance may be ob- 
tained if we consider the case of a current-carrying wire being axially strained. 

The resistance of the wire R is given by 


R= pL|A 


where L is the length of the wire, 
A = 7nr’? is the cross-sectional area of the wire, and 
p is the resistivity of the wire. 
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Taking logarithms, 
log R = logp + logL — log A 
where log A = log7 + 2logr 
Differentiating, 
d Rained A tde') dl s2idr 


But for a cylindrical wire, a longitudinal strain e, = dL/L produces a 
lateral compressive strain 
dr dL 
6,= 7S)’ 


r L 


where v is Poisson’s ratio. Therefore 
dR dp dL 
(a) = (FZ) + Fa + 
dR /dL dp /dL 
ie (qr) = 0 +294 (2/5) 
According to Bridgman (1929), for metal under hydrostatic pressure 
dp/p = c(dV/V) 


Obviously our stretched wire is not being subjected to hydrostatic pressure, 
but if we assume that Bridgman’s relationship may be applied approximately 
then we may write 


22, a NE 
5 OR Wi A 
dL 
=e7 (1 — 2r) 
Hence 
dp [aL _ 
5 Tell 2) 
dR /dL 
and R spy ICL Ae 20). ch iclL 20) 


The right-hand side of this equation gives us a measure of the sensitivity 
of a given gauge as a strain-measuring instrument and is usually replaced 
by k, the ‘k-factor’ or strain sensitivity factor for a given gauge. 

Our basic equation relating change of resistance to strain is therefore 

dR dL 


; dL 
re k Tr where strain e; = Tr 
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For most metals k = 2, so the changes in resistance are very small. With 
gauges constructed from metal wire or foil it is necessary to ‘fold’ the wire 
or foil several times so that there is a large number of closely packed parallel 


riN 
VW 


Wire Foil 
Ficure 25.1 The standard gauge. 


Galvanometer 





FicurE 25.2 Wheatstone bridge for measuring resistance changes in a strain gauge. 


strips along the axis of the gauge (Figure 25.1). This enables the gauge to have 
a reasonably high resistance. 

The small changes in resistance are measured by a Wheatstone bridge, 
illustrated in Figure 25.2. There are many modifications to this basic circuit, 
but all circuits work alike when used in the ‘null balance’ condition. In this 
case, no current passes through the galvanometer, and so 

Ri _ Re 


Ry RS 
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If R, is the gauge resistance, the bridge may be balanced by adjusting R, 
which is usually calibrated in percentage change. Null-balance bridges are 
only suitable for static measurements. For dynamic work, the bridge is 
initially balanced, and the out-of-balance current measured. 

Usually, two standard resistances (Rg and R,) are provided in the bridge, 
their resistance being the same as the strain gauges. This gives maximum 
' sensitivity. The other two resistances are both strain gauges, one stuck on the 
specimen, and the other stuck to a similar piece of material nearby. This 
provides ‘temperature compensation’, an essential factor for obtaining 
accurate results with this method. As the resistance changes due to strain are 
so small, expansion or contraction of the specimen due to slight variations 
in temperature could be read as strain. By having two gauges stuck to the 
same material at the same temperature (hence putting them close together), 
any changes in temperature will cancel out, while actual strain in the ‘active’ 
gauge (stuck to the stressed specimen) will be over and above that of the 
dummy gauge (stuck to the similar piece of material). 

The circuit may be balanced by a slide resistance wire, connecting the two 
standard gauges, and which is calibrated in percentage strain. Using this 
wire, the relative value of the standard resistances is changed, thus balancing 
the bridge. Any initial mismatch of resistance between active and dummy 
gauges may be compensated by adjusting the slide wire, or using small 
rheostats inserted in the circuit. 


Load 






Strain gauge 
locations 


Ficure 25.3 Tilted cantilever for producing combined bending and normal stresses. 


Experimental Details 


If a commercial strain gauge bridge is not available, the student may wire 
up the Wheatstone bridge as above. Usually, gauges have a resistance of 
120 ohm, so the standard resistors must have this value. The applied voltage 


6-11 
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should be 2 to 6 volt d.c. A shunt should be placed across the galvanometer 
for protection during initial adjustments. 
The gauges are supplied on a paper or thin plastic backing, with the gauge 


cauges -Y*[ >} J" 


{Cr 


<— | > 
M pds Tension M 
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Ry 
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Strain due Strain due 
to normal to bending 
load 


o% 
ad Rp 








(b) 


Ficure 25.4 Wheatstone circuits for separating bending and normal strains (a) 
circuit for bending only, (b) circuit for normal strains only. 


factor quoted by the manufacturer. One centimetre length gauges are the 
most commonly used, and are suitable for this experiment. Gauges may be 
cemented to specimens by common glues, such as Seccotine, Durofix, 
Araldite, etc. It is most important to obtain a very thin layer of cement, 
and to exclude all air bubbles. 
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(i) Using a 2:5 cm wide metal tensile specimen, cement three gauges at 
45° to each other, and another three at 60° to each other, making sure 
no gauge is aligned with the tensile axis. Cement one gauge to a similar 
piece of metal as dummy. Load the specimen and take the gauge 
readings. 

Using the formula for compound strain, and Mohr’s circle, determine 
the direction of principal strains, and their magnitudes for the 45° 
and 60° ‘rosettes’ (as the gauge combinations are known). Compare 
these results with the known values. 

(ii) The construction of the strain gauge bridge make it possible to arrange 
special circuits to separate certain types of loading, while still ensuring 
temperature compensation. Arrange for a cantilever to be loaded as in 
Figure 25.3, with strain gauges cemented on both sides of the beam. 
The bending strain alone may be measured by attaching the gauges 
to the bridge as in Figure 25.4(a), and the normal strain alone by 
the circuit in Figure 25.4(b). In this latter case, the standard resistors 
must be replaced by two dummy gauges. Load the beam, and measure 
the changes in resistance using both circuits. Measure the angle that 
the load makes with the cantilever. 

Calculate the normal and bending loads, and hence determine the 
normal and bending strains in the cantilever. Compare these with the 
measured strains. 


Discussion 


Which type of gauge rosette makes calculation of the principal strains sim- 
plest, and which is the most suitable for analysis by Mohr’s circle? Suggest a 
possible method for solving large numbers of rosette results, without doing 
the full calculation each time. Discuss the effectiveness of the dummy gauge 
for temperature compensation. 

Explain in detail how the two circuits separate normal and bending strains, 
and comment on their effectiveness, suggesting what errors may arise. 
Suggest a circuit, involving all four arms of the bridge, for measuring pure 
torsion. 
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Proving Rings 


In previous experiments we have applied the loading by means of calibrated 
weights, so that no problem existed with regard to load measurement. When 
applying loads of several tons or more weights are not very convenient, even 
when the load is magnified by a lever system. Screw and hydraulic load systems 
are used nowadays in modern testing machinery to magnify the applied load, 
which is usually measured directly by connecting a load cell behind the jaws 
of the testing machine. This direct connection avoids problems due to friction 
etc., which occur when the load measurement is made elsewhere in the testing 
circuit. ' 

Load cells are based on deflection or strain measurements using the defor- 
mation of a standard block of material, which is loaded in a simple manner 
so that there is a linear relation between load and deflection. The simplest 
device is a block of metal with strain gauges on its faces, so that change in 
length of the block due to applied load is monitored directly. 

More sophisticated instruments magnify the deflection so that less sen- 
sitive measuring devices, such as dial gauges, can be used. An obvious 
solution is to use the deflection of a beam, and a proving ring is a particular 
example of this application. These rings are widely used in load measurements, 
and are particularly useful because of their compact form and high sensitivity. 


Theory 


The deflection of a ring under diametral loading is an example of the applica- 
tion of Castigliano’s theorem, which allows deflections to be calculated from 
the strain energy. The theorem states that the deflections (or rotations) of 
the structure at the point of application of load (or moment) and in the direc- 
tion of the load equal the partial derivatives of energy with respect to load 
(or moment) at the points of application. If the strain energy can be cal- 
culated, it is fairly simple to find the derivatives. For beams, only the strain 
energy due to bending is significant, and this is given by the expression 


iz 


where the integration is taken along the length of the beam. This expression 


~~ 
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also holds for curved beams, provided that the initial radius is large com- 
pared with the dimensions of the cross-section. 

The ring is illustrated in Figure 26.1(a), with the applied load P vertical. 
The bending moment cannot immediately be obtained in terms of P and R, 
but, making use of the symmetry, let M, be the bending moment on cross- 

sections perpendicular to P (Figure 26.1(b)). There will also be a normal 
load of 4P on these cross-sections. 


Dial gauge support 


Dial gauge 


measures 
diametral 
deflection 








wl 
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(b) 


Ficure 26.1 Loading of proving ring (a) Proving ring under diametral load, 
(b) Effective loading for half the ring. 


At an angle 6 the bending moment M is given by 
M = 4PR(1 — cos 8) — Mo 


Now 


M2 
strain energy U = { Bal ds 


nr ay LB 2 
=4| {PR(1 — cos 6) ~ 2Mo}* » g 


6 
0 SEI 


From Castigliano’s theorem, the rotation of the horizontal cross-section is 
given by 


oU R 


jn 
Sa — 0) — 2M,|(—2) dé 
aM, | 2|PR(1 — cos 8) ol( ) 


=0 by symmetry. 


132 PART II STRESS ANALYSIS 


Hence 
a ; 
[ (PR — PRcos 6 — 2M,) dé = 0 
0 


PRi7 — PR — 2Mo37 = 0 
= pr llane 
sia 2 @ 
The maximum bending moment occurs when 6 = 7/2, and is 
M = PR/2 — Mo = PRi[t =0 2'@pe 
The deflection of the points under the load P is given by 
au R 
oP° 7) QET 
R2 
~ 2EI 


in 
[ 2{PR(1 — cos 6) — 2M,}(1 — cos #)R dé 
0 
ta 
[ {2PR — 4PRcos 6 + PR(1 + cos 20) — 4M, + 4M, cos 6} dé 
0 


2 
= saylPRn — 4PR + PRia — 2M, + 4My] 


wt ay “S| 
~ 4EI \ x 
We also need to ensure that the ring is never stressed past the yield point 


of the material, so that maximum stress must be calculated. This can be done 
using the basic bending equation for a curved beam, i.e. 


c M 1 | 
5-7 7F(R-R) 


and substituting the value the maximum bending moment. 


Experimental Details 


As the load-deflection equation involves the radius as a cubic term, the rings 
must be machined carefully. The material properties of the ring, e.g. Young’s 
modulus, should also be constant throughout its length. Annealing is 
advisable. 

Select ten steel rings, with diameters varying from 20 to 30cm, width 
25cm and depth about 1 cm. Five of these rings should have the same 
diameter, with depths varying from 0-5 to 20cm. The deflection of the 
rings should be measured by a dial gauge, reading to 0-0001 cm, and a 
suitable fitting for attaching this to the rings and applying the load, is 
illustrated in Figure 26.2. 

Calculate the maximum permissible applied load for each ring, using a 
safety factor of 0:5. Load each ring in turn to its maximum (a testing machine 


~ 
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of several kilonewtons capacity will be necessary) and plot the load versus 
deflection curves. Calculate the standard deviation, or otherwise estimate 





Load applied via 
ball bearings in 
conical seating 


Bolts 


Dial gauge 


FicurE 26.2 Attachment of dial gauge to proving ring. 


the linearity of this graph. Choose a load common to all the rings of varying 
diameter, and plot deflection versus cube of the radius. Repeat this process 
for rings of common radius, but plot load against the reciprocal of the depth 
cubed. 


Discussion 


Observe any significant non-linearities in the plotted graphs, and suggest 
reasons for these. Suggest any metallurgical parameters that may affect a 
ring’s performance over a period of time. Discuss the accuracy of the load 
measurements and suggest a suitably precise method for quickly converting 
deflection to load. 


Me fe | 


Structural Analysis by 
Indirect Models 





The direct determination of stresses from a model of an actual structure will 
depend on the various scaling factors involved between the model and struc- 
ture parameters. The relations between applied forces, stress, strain, dis- 
placements, etc. can be obtained from dimensional analysis, which give the 
‘laws of similitude’ between the model and structure. For a static problem 
these laws show that there are only two independent scale factors, which 
severely restricts the choice of model parameters. It is virtually impossible to 
build a model to such stringent requirements, and in practice a considerable 
variation from these conditions must be tolerated. 

If the model is to simulate the elastic behaviour of a frame structure, then 
the indirect model method is most useful. In this method, only the basic 
geometry of the model is scaled to the structure, and the loading is applied 
by arbitrary deformations. The deflections of the model are proportional to 
its influence lines, from which the bending moments, shear forces, etc. and 
hence stresses may be calculated. For complicated redundant structures, this 
method has been used as an effective design tool, but now it can be replaced by 
numerical procedures of solution using digital computers. 


Indirect Model Analysis 


Influence lines are curves that show the variation of reaction, bending — 
moment or shear force at a point in a structure as a load moves across it. 
For a given load, a structure will possess any number of influence lines, 
depending on the section under consideration. 

The experimental determination of influence lines is based on the Maxwell 
reciprocal theorem and the Muller-Breslau theorem. Maxwell’s theorem states 
that the deflection of an elastic structure at point A due to an applied load at 
point B equals the deflection of the beam at point B when the same load is 
moved to point A. The directions of load and displacement must correspond 
for this theorem to hold, but in most applications the directions are either 
vertical or horizontal. 

Muller-Breslau postulates that in a statically indeterminate structure the 
deflection of the centre line of the structural components, produced by 


~ 


STRUCTURAL ANALYSIS BY INDIRECT MODELS 135 


displacing one of the supports in the direction of a redundant reaction, 
defines the influence line for that reaction to some scale. This theorem also 
applies to rotations and shearing displacements (produced by cutting the 
frame) made on the frame. 

The application of these theorems is best explained by reference to a 
simple indeterminate structure. Consider the propped cantilever in Figure 
-27.1(a). We wish to determine the influence line for the support reaction R, 


Unit load 


(a) 









a 82 
Pom load 


(b) 


Unit load 


ee ee por 


(c) 


Ficure 27.1 Application of indirect model theorems to a propped cantilever, 
(a) Propped cantilever with unit load in position, (b) Support removed and 
replaced by unit load, (c) Deflection of beam due to the unit load in its original 
position. 


i.e. we wish to find the curve showing the variation of the reaction R when a 
load (taken as unity in value) is moved across the cantilever. The support is 
removed, and replaced by the unit load in Figure 27.1(b) and this produces a 
displacement 6, and a corresponding displacement 6, at the original position 
of the unit load, which can be taken as anywhere along the beam. From 
Maxwell’s reciprocal theorem, the deflection 6, equals the deflection at the end 
of the beam (with the support removed) when the unit load is replaced at its 
original position, as in Figure 27.1(c). But this deflection is zero when the 
support is in position, due to the applied load R, and as a unit load produced 
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deflection 6, at this position, the reaction would produce a deflection Réz. 
Thus the two opposing deflections cancel out, giving 


Ré, = Oy 
R = 6,/62 i: 


By measuring the deflections of the beam along its length, and dividing by 
the displacement due to the unit load applied to the end of the beam, we 
obtain the influence line for R. Experimentally this procedure is very simple, 
and any number of supports can be introduced without increasing the labour. 

The same principle can be applied to the determination of bending moments 
and shear forces, but this is more difficult in practice as arbitrary rotations and 






Balance beam 


Model member 


FiGurE 27.2 Bending moment deformeter. 


shear displacements must be introduced into the structure. Special defor- 
meters have been introduced for this purpose, for use with models using thin 
plastic strips which make it unnecessary to cut the structure. A bending 
moment deformeter is illustrated in Figure 27.2. 

The method is quite accurate providing (a) the structure is linearly elastic, 
(b) there are no stability effects, and (c) there is no gross change in geometry 
under load. This latter point infers that all the displacements are small, 
which experimentally is not convenient. For most structures, it is possible to 
use displacements that can be measured easily with a ruler without introduc- 
ing significant errors. Plastics are the most convenient model materials 
because of their low modulus, but they also creep considerably at room tem- 
perature. For accurate work this can be overcome by applying the displace- 
ments through creep compensators, which are made of the same material as 
the model. 


Experimental Details 


The propped cantilever is an ideal subject for study, as it is basically simple 
and several supports can be added quite easily. The beam can be a Perspex 
strip, at least 30 cm long, about 3 mm thick, and 1 cm wide, or a thin metal 
strip about 0-1 mm thick may be used if available. The strip must be clamped 
firmly between two metal jaws, which are mounted on a pivot and slotted 
into a drawing board so that the edge of the beam is against the board, 
as in Figure 27.3. 
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The first set of experiments will determine the influence lines for the 
support reactions as already described. 

Fix the pivot rigidly in its mount, and first determine the influence line for 
the reaction of a single support placed at the beam’s free end. The arbitrary 
displacement can be made by knocking panel pins into the board to support 
the beam. The deflected shape is best drawn on graph paper, which is pinned 
on the board prior to loading. The centre line of the beam should really 
be drawn, but as the beam is thin, its bottom or top may be used. It is also 
more accurate to deflect the beam by equal amounts above and below the 
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FicurE 27.3 Arrangement of the cantilever on a drawing board. 


neutral position. The displacements are then measured between the two 
deflected positions. 

Next add a further support to the centre of the beam, and measure the 
influence lines for the reactions in each support whilst the other is in position. 
The supports may be fixed by using panel pins on either side of the beam, 
but make sure the fixing is firm. 

Finally, determine the influence line for the bending moment at the encastre 
end of the beam with only a single support at the other end. In this case, the 
pivot clamp must be released and the clamp given an arbitrary rotation. This 
may be measured by fixing a long pointer to the clamp. The bending moment 
at the clamp due to a unit load moving along the beam is given by the ratio 
of beam deflection divided by clamp rotation (measured in radians). 

Plot the influence lines for each of the three cases tested. The reaction and 
bending moment for the single support cases can be calculated fairly simply, 
as the deflection at this point must always be zero. Obtain the expression for 
the deflection of the free end due to unit load a distance x from the clamped 
end, and equate this to the expression for deflection due to a load R supporting 
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the end of the beam. This will yield R in terms of the beam constants (length, 
modulus, moment of inertia) and the distance x, and the influence line 
for Rcan be calculated. Use this expression to determine the bending moment 
influence line, and compare both calculated results with experiment, by 
plotting about six calculated results on the experimental graphs. 


Discussion 


Discuss what sources of error may be significant in your experimental results, 
and why it is more accurate to measure displacements between beams 
deflected in opposite directions. Suggest how you would use influence lines 
to determine maximum stresses in your structure, and hence aid the design 
process. 


28 


The Photoelastic Method of 
Stress Analysis 


Many transparent materials that in an unstressed state are not naturally 
double-refracting become so under the influence of stress. Such an effect 
is known as artificial or stress birefringence, and the magnitude of the optical 
effect produced by stress is usually a direct linear function of the stress 
magnitude. 

The first recorded investigation of this effect was made by Sir David 
Brewster, who discovered that a simply loaded plate of glass exhibited the 
double refraction effects normally associated with a uniaxial crystal, the axis 
of load being equivalent to the crystal’s optic axis. 

The theoretical treatment of this phenomenon was developed indepen- 
dently by both Neumann and Maxwell, although Neumann related the 
optical effects to strain, whereas Maxwell expressed his equations in terms of 
stress. Whether, in fact, the observed phenomena are strain-based or stress- 
based only becomes of importance, however, when problems of a non- 
elastic nature are considered. 

The use of photoelastic models for the stress analysis of structures by 
analogy was first developed into a practical technique by Coker and Filon, 
who in over twenty years of research at University College, London, laid the 
foundation of modern methods of photoelastic stress analysis. Their classic 
treatise on the work done during this period is still the most authoritative 
publication in this field. 


The Optical Basis of the Photoelastic Effect 


It is convenient for our purposes to consider light as a transverse electro- 
magnetic vibration. Visible light has a wavelength between approximately 
4000 and 8000 A (1 A = 10-® cm). 

Normal light is unpolarized in that there is no predominant plane of 
vibration for the electric vector of the wave. Light can be polarized however, 
by successive reflections or refractions at critical angles, by electrostatic 
or electromagnetic fields, or by passing the light through various types of 
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crystal. Modern practice, however, is to use ‘Polaroid’ sheet, a very efficient 
polarizing medium, which can be considered as acting on randomly vibrating 
light waves in the manner shown schematically in Figure 28.1. When light 
polarized in this way passes through a transparent substance that exhibits 
stress birefringence, it is split into two orthogonally polarized components, 
each parallel to a direction of principal stress. These components propagate 
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FiGurE 28.1 The photoelastic effect—plane and circularly polarized light. 


through the material at different velocities. For light of a fixed wavelength 
(monochromatic light), the velocity of each of the components is propor- 


tional to the magnitude of the principal stress lying in its plane of vibration, 
jf 


Vi, £04; Vo Oo 


Thus, on emerging from the photoelastic medium, the two components will 
have suffered a relative change of phase or ‘relative retardation’ that will be 
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proportional to (V; — V2). But as (V,; — V2) is in turn proportional to 
(0, — o,) we can express the relative retardation 6 as 


6 = Ct(o, — oe) (28.1) 


where C is the stress-optic coefficient and ¢ is the thickness of the specimen. 

If we view the emergent light through another sheet of Polaroid crossed 
with the first and called an ‘Analyser’ (Figure 28.1), we see only the horizontal 
components of these two rays. When the relative retardation is such that 
the two components are vibrating exactly in phase with one another, mutual 
extinction will occur, producing a black fringe for every point in the photo- 
elastic structure which has a similar value of principal stress difference 
(o, — oz), assuming that C and ¢ in equation (28.1) are constant over the 
structure. 

A full mathematical analysis of the plane polariscope system shown in 
Figure 28.1 shows that, for light of a given wavelength A, the transmitted 
intensity of light emerging from the analyser is proportional to 


t 
sin? 2a sin? “ ay (28.2) 


where a is the angle between direction of vibration of polarized light and 
principal stress o;, ¢ is the thickness of photoelastic medium, A is the wave- 
length of monochromatic light, and n, and nz are the refractive indices of the 
two components vibrating in the directions of 0, and og respectively 
(n, = V/V,, ng = V/V2, where V is the velocity of light in air). 

Thus we can see that there are two conditions for extinction (zero intensity) 
of the emergent light. 


(i) Everywhere on the photoelastic structure where a = 0° or 90° (ie. 
either o, or o2 aligned with polarizer axis), extinction will be produced 
irrespective of the particular wavelength of light being used. Thus a 
black fringe will be seen across the structure wherever the above 
conditions apply. This fringe is called an isoclinic, and it gives us the 
locus of points in the structure at which the principal stresses have a 
given direction (that of the polarizer and analyser axes). Hence by 
rotating the structure relative to the polarizer and analyser and so 
bringing the isoclinic fringe to a given point, it is possible to determine 
the directions of principal stresses at every point on the structure. 

(ii) Extinction of emergent light will also occur when 

(wt/A)\(n, —1m2)=0 or iwmi=1,2,3,..., 
But t(ny paps Ng) = 6 
where 0 is the optical path difference or relative retardation of one com- 
ponent on the other. Therefore the condition for extinction is given by 
6 
c =0 of iz (28.3) 


or 6=0 or td 
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Thus when the difference in optical paths of the two components is 
zero or an integral number of wavelengths, extinction will occur for a given 
wavelength 1. Note that 


6 = tm, — na) 


ag li 7) 
pad We eae 


Ve Sage “| 
pis ( VV, 


Vt 
=P (V2 — Vi) 


where V, is the velocity of the wave in the unstressed medium, V is the velocity 
of the wave in air, and V, and Vz are the velocities of the wave in the planes of 
o, and o, as previously defined. But in equation (28.1) we have already 
defined relative retardation 6 as equal to Ct(o, — o,)so wesee that as (co, — o2) 
and hence (7, — m2) increases in magnitude from zero there will be successive 
extinctions at equal intervals, producing fringes of zero, first, second, third 
order, and so on. These fringes are called isochromatics and provide a means 
whereby the principal stress difference or maximum shear stress $(o, — o2) 
may be determined. 

If, instead of using monochromatic light, we use white light, we see that as 
isoclinics are independent of wavelength, they will still be seen as black 
fringes. Isochromatics, however, being dependent on wavelength (equation 
28.3), will appear coloured in white light, the colour being given by the 
complementary colour observed in white light when that particular wave- 
length is removed by the process described in (28.2). 

If we observe, in white light, a linearly increasing stress-gradient such as 
can be obtained with a cantilever, the colour sequence, starting with black 
(which we may term a ‘zero-order’ isochromatic), will be as shown in 
Figure 28.2. It will be noted that as we obtain higher order isochromatics 
their tint changes, as there is more than one value of i that will give extinction 
in a higher order spectrum. For instance, if 6 is 30000 A, A(= 6/i) can 
equal 7500, 6000, 5000, 4280, and 3750, producing five separate extinctions 
in the visible spectrum. 

In very high orders (7 or 8) the spectrum will be interspersed at even 
intervals with extinction bands from one end to the other, the net effect to 
the human eye being indistinguishable from normal white light. Thus as 
fringe orders increase, the fringes become more and more pale and eventually 
‘wash out’ completely if white light is used. If monochromatic light is used, of 
course, no such barrier exists to the observation of high order fringes. 

To sum up, in plane-polarized white light we observe both isoclinics 
(black fringes) and isochromatic (coloured fringes). 

If the photoelastic structure is rotated in relation to polarizer and analyser, 
the isoclinics will move (relative to the structure) but the isochromatics will 
not. The movement of the isoclinics across the structure is used as a means of 
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FiurE 28.2 The use of white light in photoelasticity. 
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determining the direction of principal stresses at a point. The isochromatics 
are used to determine the magnitude of the principal stress difference or 
maximum shear stress at a point, as each colour observed corresponds to a 
certain stress magnitude (Figure 28.2). 

Accurate measurements of stress magnitudes are made by means of various 
compensation techniques, either goniometric (Tardy or Senarmont) or by 
using optical compensators. 

In many cases it is found more convenient to use circularly polarized light 
when determining stress magnitudes. A circularly polarized polariscope 
arrangement is shown in Figure 28.1. 

Two quarter-wave plates are inserted as shown with their fast axes crossed 
and at 45° to the polarizer and analyser axes. The first 2/4 plate splits the 
plane-polarized light into two orthogonal components of equal magnitude, 
out of phase by a quarter of a wave-length. The resulting electric vector 
produced by these two components thus rotates in a circle as shown sche- 
matically in Figure 28.1. On passing through the structure, the light is in 
general elliptically polarized as shown, and is then finally resolved by the 
second 4/4 plate and analyser. 

A full mathematical analysis of the system shows that the intensity of light 
emerging from such a system is zero when 


t(ny — Nz) = 0 or t,t ie, Be 


This is the condition for isochromatics and thus in the circular polar- 
iscope the isochromatic pattern can be seen without the superposition of 
isoclinics. 

In order to return to plane-polarized light, the 4/4 plates are simply removed 
from the system. It should be noted that when white light is used in conjunc- 
tion with the circular polariscope, then because the quarter-wave plates 
correspond to 4/4 for only one specific wavelength, the rest of the spectrum 
on either side of this wave-length will be elliptically polarized to a greater 
or lesser extent. As long as quantitative readings are taken relative to the 
colour corresponding to this particular wavelength, the readings will still 
be valid. 


Applications 


The photoelastic method will be used to study a tensile specimen, a beam 
under four point loading, and a tensile specimen with a hole. 


Experimental Details 


If a polariscope is not available, a simple instrument may be constructed from 
two pieces of Polaroid, two quarter-wave plates (also available from the 
Polaroid Corporation), a one metre V-type optical bench with four saddles, 
and a simple loading frame that can sit astride the bench, as shown in Figure 
28.3. Two 7-5.cm diameter 15 cm focus plano-convex lenses, a projection 
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lens (an enlarging lens will do for this), a screen, and a 500 W projector 
lamp in its housing complete the equipment. 

The Polaroid should be cemented to the two plano-convex lenses, which 
should be mounted in holders capable of rotating about the optic axis, and 
which are calibrated in degrees, the zero position corresponding to the polariz- 
ing axis. The lenses should be placed on either side of the loading frame, with 
the Polaroid facing the frame. The projector lamp is placed at the focus of 
one lens to produce collimated light, and the projector lens and screen are 
placed on the same side of the other plano-convex lens, their positions being 
adjusted to give a sharp image of the model in a frame of suitable size. 
The quarter-wave plates are inserted with their fast axis at 45° to the Polaroid 
axis (and between the Polaroid and model) and mounted on separate saddles, 
so that they can be easily removed. 





Ficure 28.3 Optical system of polariscope: S is the projection lamp, L are plano- 
convex lenses, P is a Polaroid, Q are quarter-wave plates, M is a model, R the 
projection lens and C the screen. 


The models may be constructed from epoxy resin or polyurethane. Cast 
sheets of epoxy resin may be obtained commercially. The models, which should 
be 6 mm thick, may be cut to the relative dimensions shown in Figure 28.4. 
All cutting should be done on a band-saw with a carefully filed finish to 
prevent residual stresses due to overheating. Similarly, all holes should be 
drilled under size and then reamed out. 

The specimens should be about 6 mm thick and 13 mm wide (allowing 
extra width at the ends of the tensile specimen for attachment to loading frame). 
These dimensions are important, to avoid torsional effects and to obtain a 
reasonable number of fringes. 

First examine the tensile specimen under increasing load, observing the 
isochromatics only. Note the relation between colour and load, and so 
obtain the ‘fringe factor’ (load necessary to produce one fringe) for the 
material you are using. In white light, the position of a fringe may be taken 
as the region between the red and blue parts of the spectrum, known as the 
‘tint-of-passage’. 

Next observe the isochromatic fringe pattern in the region of the beam 
undergoing pure bending, and sketch this pattern when there is sufficient 
load to produce at least three (and preferably more) fringes in the field of view. 
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Remove the quarter-wave plates, and rotate the polarizer and analyser 
together to observe the positions of the black isoclinics, again restricting your 
observations to the region of pure bending only. 

Finally observe both isochromatic and isoclinic patterns for the tensile 
specimen with a hole drilled in its centre. Sketch the isochromatic pattern 
around the hole when the load is sufficient to produce at least three fringes, and 
also the isoclinic pattern when the polarizer axis is parallel to the tensile axis. 
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FIGURE 28.4 Dimensions of photoelastic models. 


Compare the fringe patterns obtained from the tensile specimen and beam 
with the patterns you would expect from the stress distribution in each of 
them (remember that isochromatics are contours of maximum shear stress, 
and that isoclinics represent points where the principal stress directions 
are parallel to the polarizer and analyser directions). 

Examine your sketches of the ‘hole’ specimen, and locate the position 
of largest maximum shear stress. Now examine the isoclinic pattern, and 
determine the principal stress directions at this point (these directions can 
also be inferred from the specimen geometry—how ?) From this information 
determine the values of the principal stresses, and hence the stress concentra- 
tion factor (i.e. the ratio of maximum principal stress around the hole to the 
mean tensile stress). 


Discussion 


Discuss the usefulness of the photoelastic method for experimental stress 
analysis. Suggest one possible method for determining the values of each 
principal stress from photoelastic data. 
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Measuring Deflections by the 
Moiré Effect 





The Moiré effect is the name given to the fringe patterns observed when two 
straight-line grids are placed in contact, but not accurately aligned. The effect 
is very common in everyday objects, such as loosely woven cloth, sets of 
railings, television, etc. 

The phenomenon has been known for a long time, but it is only in the last 
two decades that any attempt has been made to use it in technology. In 
experimental stress analysis it has been utilized to measure surface strains, 
flexure and deflection of plates. This experiment, which requires the minimum 
amount of equipment, will demonstrate the latter. 

The effect is illustrated in Figure 29.1(a), where two grids of slightly 
different pitch are aligned with their grid lines parallel. This produces fringes 
parallel to the lines where the lines of one grid obstruct the spaces of the 
other grid. Thus the fringes represent the difference of one pitch between each 
grid, i.e. a fringe occurs when one grid ‘catches up’ on the other by one line. 
Measurement of the fringe spacing gives the difference between the two 
pitches. A similar effect occurs when two grids with identical pitches are 
misaligned as shown in Figure 29.1(b). A more general Moiré pattern is 
made up of a combination of (a) and (b). 


Theory 


Deflections are measured by projecting a grid onto a diffuse reflecting sur- 
face, and then observing the shadow of the grid through a similar grid. This is 
achieved quite simply by placing a grid in close contact with the specimen 
and illuminating it with light at an angle i to the normal, thus projecting the 
grid onto the surface. If viewed normally, then any variations in the depth 
of the specimen beneath the grid will produce a fringe pattern, as illustrated 
in Figure 29.2. From this diagram, we see that the fringes represent contours 
of constant displacement, where the contour interval is 


h = pitani 
where p is the grid pitch. 
Using a grid of 100 lines/cm and an illuminating angle of 20°, the contour 
interval is 0-025 cm which is a useful range for beam and plate problems. 
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Ficure 29.2 Principle of the Moiré ‘shadow’ method: each fringe position corres- 
ponds to an increase of one pitch (p) between the light-source ray and the camera 
ray, so that the fringes are contours of constant depth h from the grid. 
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Ficure 29.3 Anticlastic effects (a) anticlastic curvature in thick beam caused by 
bending, (b) bending of a thin, wide plate. 
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greater bending stresses to counteract the Poisson effect, and this produces an 
apparent increase in the beam stiffness of 1/(1 — »?) if simple beam theory is 
used. 


Experimental Details 


Use six Perspex rectangular-section beams, about 20 cm wide with depths 
ranging from 5 cm to 1 mm. Paint the top surface with a matt white finish, 
and load the beams to obtain a section under pure bending. The surface of 
this section should be flat to within 5 x 10-* cm, and if necessary these 
surfaces should be machined beforehand. 

A Ronchi grating, with between 80-120 lines/cm should be placed on top 
of the beam, and illuminated with a Photoflood bulb at least 60 cm away, 
and about 30° to the specimen normal. The surface should be observed 
normally, either with a camera, or a cathetometer with which measurements 
can be made directly. It is not necessary to shut off any other lights, providing 
they are a considerable distance from the experiment. The lighting should be 
adjusted to get the best possible fringes, but try to keep the angle as large as 

ossible. 
: As the angle of the illumination cannot be determined directly, an initial 
calibration test must be done using a beam of similar dimensions to the test 
pieces but machined with a 200 cm radius of curvature. 

If this is placed in position with the grid on top, a reasonable number of 
fringes will appear and the illumination may be adjusted for maximum 
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FicurE 29.4 Geometry of a circular beam for calibrating Moiré-apparatus. 
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fringe contrast. The fringe separation can be measured and as the deflections 
lie on a circular arc, Figure 29.4, the depth at any point along the beam can 
be calculated from the formula 


2Ry = x? 
where x is the distance along beam from uppermost part of arc, but 


= Nhe N= 
NE = N ean 


where N is the fringe order, or integer. 


Thus by plotting N against x?, we get a straight line graph of slope tan i/2Rp, 
so that tan i can be calculated. 

The deepest beams will have to be loaded up to maximum bending stress 
to obtain any fringes at all. Sketch the fringe pattern for each beam, and 
measure the fringe positions along the centre and the two edges of each beam. 
The shallow beams should be loaded to obtain a reasonable number of fringes 
not more than twenty. 

Remembering that pure bending produces a circular deflection, use the 
measurements along the beam to calculate R,, and hence the value of Young’s 
modulus for each beam, from the simple bending equation 


E|R, = M,|I 


For the deepest beams only, calculate the transverse curvature using the 
three measurements made on any one fringe and assuming that this trans- 
verse curvature is circular. Compare the two radii of curvature to find a 
value of Poisson’s ratio. 


Discussion 


Comment on the variation of Young’s modulus that you obtain, and also on 
the Poisson ratio value. Do your stiffness values increase by as much as 
1/(1 — »?)? 

What further uses do you foresee for this application of the Moiré effect ? 
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Brittle Lacquer Methods 


If a structure is coated with certain types of lacquer which harden into a 
very brittle thin coating, then when the structure is loaded cracks will appear 
in the coating in a direction perpendicular to the greater principal tensile 
stress. The distribution of cracks will indicate graphically the directions 
of the surface stresses and will show up areas of stress concentration. At 
this stage it is common practice to use the information obtained as a means of 
determining where to put strain gauges, which are then used to determine 
the actual strain magnitudes. However, provided that the lacquer is applied 
under closely controlled conditions, it is possible to determine approximate 
values of principal stress magnitude from the load level at which individual 
cracks are initiated. 

There are a number of materials which can be classed as brittle lacquers 
and, in general, they possess the following properties: 


(a) they develop well-defined cracks at strains below the yielding strain 
of the structure; 

(b) the cracks, once developed, remain open when the straining load is 
removed; 

(c) the initial cracking response to strain is independent of the thickness 
of the coating within reasonable limits; 

(d) the density of the cracks at a given strain varies with the thickness of 
the coating—the thicker the coating, the fewer the number of cracks 
per centimetre for a given strain; 

(e) the coating does not flake at reasonable strain levels. 


A wide variety of lacquers are suitable for this work, and even whitewash 
is used in some engineering works to detect any large plastic strains in struc- 
tures. Unfortunately, most lacquers are very sensitive to changes in tempera- 
ture and/or humidity, and this makes their quantitative use extremely difficult. 
They detect only tensile strains, but this problem may be overcome by pre- 
loading compressed specimens prior to the lacquer application, hence detecting 
the apparent tensile strains on unloading. 

Despite these difficulties, the brittle lacquer technique is a widely used tool 
in machine stress analysis and, combined with resistance strain gauges, it 
can solve many practical stress analysis problems, particularly on curved 
surfaces which are difficult to analyse. 
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Application 


A brittle lacquer will be used to examine the effect of stress raisers, such as 
holes and notches, in Perspex models. Any model may be used, but useful 
examples are provided by stress distributions near a hole and a notch in 
tensile specimens. An attempt should be made to obtain quantitative 
measurements by spraying a Perspex cantilever with lacquer for use as a 
calibration bar. 


Experimental Details 


A useful lacquer is made up as follows: 


160 g LiOH, H,O 
480 g H,BO, : 
454 ¢ distilled water to make about 800 ml solution. 


0-7 g Na fluoresceinate 


The result is lithium tetraborate, as developed by Hickson at the Royal 
Aircraft Establishment who states ‘Ingredients should be the purest obtain- 
able, and all vessels should be very clean and free from crystallization nuclei. 

The lithium hydroxide is added to the water at 60 °C and stirred until 
most has gone into solution. Add the boric acid while stirring vigorously, 
continuing the agitation until complete dissolution has been effected. 
Dissolve the dye. Filter under vacuum through a No. 4 sintered glass 
filter.’ 


This lacquer should be applied at constant temperature and humidity, as 
it is very sensitive to ambient conditions. It may be sprayed or painted, and is 
non-toxic. It dries rapidly at room temperature. 

Coat the specimens and calibration bar at the same time, ensuring the same 
even coating in each case. Allow to dry for one hour. 

Test the coating by applying increasing load to the cantilever bar. The 
strain at which the first crack appears is the ‘threshold strain’, and should be 
between 5 and 10 x 10-%. However, as the calibration bar is the same 
material as the specimens, it is more useful to calibrate in terms of stress. 
As the load is increased, more cracks appear perpendicular to the maximum 
principal stress direction, and will be most closely spaced in the region of 
highest stress, i.e. the stress raiser. Thus by counting the number of cracks/cm 
over a small length, perpendicular to the cracks, an estimate of the maximum 
principal stress can be made. 

However, the most useful method of quantitative analysis is to increase 
the load in equal increments, and mark the position of the threshold cracks 
for each increment (Figure 30.1). In this way, a pattern of contours of equal 
tensile stress intervals can be built up, the stress level being highest at the first 
contour to appear. 

Increase the load in equal increments for both calibration bar and 
specimens, and plot the contours as described. 
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Compare your stress maps of the calibration bar with the calculated stresses 
by plotting the two graphs of stress against length. 


oak 





== —=— Stress 
contours 


FicurE 30.1 Stress contours from brittle lacquer cracks, by the method of equal 
load increments, 


From the specimens, deduce the position of maximum stress, and estimate 
the stress-raising factor. 

Remember that the brittle lacquer technique is not very reliable, and you 
will be lucky to be within 20% of the theoretical solution. 


Discussion 


Discuss your results, and suggest possible faults in your particular technique. 
Comment on the usefulness of the method. 
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